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We consider two-dimensional Hamiltonians on a torus with finite range, finite strength interactions 
and a unique ground state with a non-vanishing spectral gap, and a conserved local charge, as 
defined precisely in the text. Using the local charge operators, we introduce a boundary magnetic 
flux in the horizontal and vertical direction and evolve the ground state quasi-adiabatically around 
a square of size one magnetic flux, in flux space. At the end of the evolution we obtain a trivial 
Berry phase, which we compare, via a method reminiscent of Stokes' Theorem, to the Berry phase 
obtained from an evolution around a small loop near the origin. As a result, we prove, without any 
averaging assumption, that the Hall conductance for interacting electron systems is quantized in 
integer multiples of e^ /h up to small corrections bounded by a function that decays as a stretched 
exponential in the linear size L. Finally, we discuss extensions to the fractional case under an 
additional topological order assumption to describe the multiple degenerate ground states. 



At low temperature, the Hall conductance of a quantum system can be quantized to remarkable precision. While this 
is an experimental fact, and the essential ingredients of our intuitive understanding of it were provided by Laughlin[5], 
there is still no fully satisfactory mathematical proof of why this happens. The main approaches theoretically are 
either the Chern number approach jH [5], which relies on an additional averaging assumption, and the noncommutative 
geometry approach[6j which is only applicable to non-interacting electrons. In this paper, we present an approach 
which avoids averaging and is applicable to interacting electrons, and can even be extended to handle fractional 
quantization. 

Before we formulate our main result, we discuss our setup of the Quantum Hall Effect. We will consider a discrete, 
tight-binding model of interacting electrons with orbitals centered at sites on a torus T. In particular, T is obtained 
by joining the boundaries of a finite [I, L] x [l,L] subset of Z^. At each site s e T, we introduce the charge operator 
Qs with eigenvalues 0, 1, . . . , gmax representing the states with the respective charge occupying the site. Let Az be 
the algebra of observables associated with Z <zT. From a general point of view, we are interested in properties of the 
ground state of the Hamiltonian Hq = X]zct*^('^)' where the following assumptions are satisfied: 

1. The interactions ^{Z) have finite strength and finite range. Formally, 3 J, i? > with L > 2R such that: 

(a) $(Z) = $(Z)t e Az, 

(b) sup,^tEz3s\\HZ)\\<J, 

(c) VZ C T, diam(Z) > R =^ $(Z) = 0, 

where 



diam(Z) = max {d(si,S2)}, d{si, S2) = \x{si) - x{s2)\ (mod L) -|- |y(si) - y(s2)| (modi), 

with x(s), y{s) denoting the x and y coordinates of site s, respectively. 

2. The Hamiltonian Hq has a unique, normalized ground state, which we will denote by I'I'o), and a spectral gap 
7 > to the first excited state. We use Pq — |^'o)(^o| to denote the projector onto the ground state. 

3. Finally, the total charge Q = X^sgt conserved, so that [Q,Ho] = 0. 

Our main result is the following: 

Theorem 1. Let Hq be a Hamiltonian satisfying the above assumptions. Then, the Hall conductance axy of |^o); 
defined in (127) for a system of interacting particles described by Hq, satisfies the quantization condition 
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for some n € N, for all L such that the following inequality is satisfied: 



(2) 



In e'^/h denotes the square of the electron charge divided by Planck's constant. The quantities C,C' are numeric 
constants of order unity, and v = 132e (1 + R)^J < 360(1 + R)^J. 

This implies that, for fixed J,R,j and g,nax, the Hall conductance is exactly quantized to an integer multiple of 
e'^/h in the thermodynamic (L oo) limit. Note that, for fixed R, J,'y,qmax, the inequality ^ is satisfied for all 
sufficiently large L. 

Before giving the proof, we discuss the applicability of our assumptions to physical experimental systems. The 
first assumption includes only finite range interactions. While there are long range Coulomb interactions in real 
experimental systems, the screening of Coulomb interactions may justify this assumption. Moreover, the case in 
which the torus is not exactly square, but has an aspect ratio of order unity, can be handled by combining several 
sites in one direction into a single site to make the aspect ratio unity; in fact, with minor changes in the proof, a 
polynomial aspect ratio can be handled at the cost of polynomial changes in the prefactor of ([T]). We omit this case 
for simplicity. Our techniques can also be extended to the case of exponentially decaying interactions, which we also 
omit. 

We note here that sites s, s' € T with y(s) = y{s') and x(s) = L, x(s') = 1 are neighboring sites on T. The same 
is true for s,s' € T with x{s) = x{s') and y{s) = L, y{s') = 1. Thus, interactions ^{X) in Hq on subsets X which 
cross the boundaries x = and y = are not necessarily vanishing. Finally, we note that without loss of generality, 
we may assume that each term ^{Z) of the initial Hamiltonian commutes with the total charge Q, by substituting 
4>(Z) with the averaged interaction which has the same range R and strength still bounded by J: 



Our proof is closely related to the Chern number approach jH [5] , so we review it and then contrast our technique. 
In the Chern number approach, one computes the Hall conductance of the above system of interacting electrons by 
introducing magnetic fluxes through twists 9^ , 9y at the boundary of the torus. A time-dependent flux 9^ creates an 
electric field in the x direction and the flux 9y measures the current in the y direction, and thus the Hall conductance 
can be identified with a certain current-current correlation evaluated at the ground state IvPo)- In the Chern number 
approach, one must assume that there is a unique ground state for all 9x,9y. If the ground state is adiabatically 
transported around an infinitesimal loop near 9^ = 9y = 0, the state acquires a phase proportional to the area of the 
loop multiplied by the Hall conductance. That is, the Hall conductance is equal to the curvature of the connection 
given by parallel transport of the ground state. The average of the Hall conductance over the torus is then equal to an 
integer, the Chern number. Unfortunately, this argument requires a non-vanishing gap for all 9^, 9y and only provides 
a quantization result for the averaged Hall conductance and does not provide a quantization result for 9x = 9y = 0. 

Our proof avoids both of these assumptions, while only requiring the presence of a gap at 9x ~ 9y = 0. Roughly, 
our argument proceeds as follows. Given any state and any path in parameter space, we define the quasi- adiahatic 
evolution, introduced in [T], of that state along the path. Various different quasi-adiabatic evolutions have been defined 
in [HlllIT^]; the evolution we choose here is yet a different one. The essential features of quasi-adiabatic evolution are 
that the evolution of a state be described by an approximately local operator acting on that state and so long as the 
path is in a region of the parameter space with a sufficiently large spectral gap, the ground state of the Hamiltonian 
at the start of the path is quasi-adiabatically evolved into a state close to the ground state of the Hamiltonian at the 
end of the path, up to multiplication by a phase. 

Our proof then rests on four results. First, in corollary |3j we show that the quasi-adiabatic evolution of the ground 
state around a small (but non-infinitesimal) square loop near the origin in parameter space (9x = 9y = 0) gives a state 
which is close to the ground state multiplied by a phase, which we relate to the Hall conductance at 9x = 9y ^ 0. 

Second, we show that quasi-adiabatic evolution of the ground state following the path starting at (0,0), going to 
{0,9y), then to {9x,9y), evolving around a small loop at the given {9x,9y), and then returning to (0,0) again gives 
a state close to the ground state up to a phase. We do not require any lower bound on the gap on the path in this 
step, and if the gap does become small then the quasi-adiabatic evolution around this path may not approximate the 




Sketch of the main argument. 
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adiabatic evolution. Rather, the proof that we return to the ground state depends on locality estimates. Crucially, 
we show in lemma [6] that the phase produced by evolving around this path is approximately the same for any 9x,0y 
Third, we show, in a procedure reminiscent of Stokes' theorem, that the quasi-adiabatic evolution around a large loop 
in parameter space, 

A : (0,0) -> (2^,0) -> (2^,27r) -> (0,27r) (0,0), 

can be exactly decomposed into the product of a large number of quasi-adiabatic evolutions around small loops as is 
described at the beginning of this paragraph (see also. Fig. |2|. Fourth, we show in corollary |4] that quasi-adiabatic 
evolution of the ground state around the large loop returns to a state which is close to the ground state with a trivial 
overall phase. 

Combining these four steps, the trivial phase produced by quasi-adiabatic evolution around the large loop is close 
to the product of the phases around each small loop. Since the phase produced by evolution around each small loop 
is approximately the same, we can then show that the phase around a small loop near the origin, raised to a power 
equal to the number of small loops, is close to imity. Since this phase is related to the Hall conductance, this gives 
the desired proof of quantization. 

The most difficult technical step is the second one. The first step relies only on the fact that quasi-adiabatic 
continuation for systems with a large gap approximates adiabatic evolution. The fourth step is similar to the proof 
of the Lieb-Schultz-Mattis theorem in higher dimensions [T] . We use energy estimates and additional virtual fluxes 
(j>x = —dx; 4>y = —0y introduced at x = L/2 + 1 and y = L/2 -I- 1 to show that the quasi-adiabatic evolution of the 
ground state along any of the four sides of the square (for example, from 9^ = 0y = to — 27t, 6y = 0) gives a state 
which is close to the ground state up to a phase. Then, we show that these phases cancel between evolutions along 
opposite sides of the square. 

Lieb- Robinson bounds play a critical role in our proof. These bounds were first introduced in and extended 
in [9l [To] . Nachtergaele and Sims, in [10] , gave the important extension to general lattices including exponentially 
decaying interactions. The most recent, tightest, and most general estimates, which we use, are from [7]. The need 
to use Lieb- Robinson bounds does currently limit us to lattice Hamiltonians; the extension to fermions moving in 
would require using unbounded interactions and an infinite dimensional Hilbert space and is currently beyond our 
techniques, although work such as [T7] is an important step towards this. 

We use C throughout the text to refer to various numeric constants of order unity. We use the phrase "for sufficiently 
large L" to indicate that L is large enough that Eq. ^ is satisfied. There are three reasons that we have a lower 



bound on L. The first reason is that in many places, such as Eq. ( 120 1, we construct subsets of the lattice whose size 
is a fraction of L; this first reason in fact only requires that L be at least a constant times R. The second reason is 
that many of the error estimates are a sum of error terms which have different behaviors in L, and we want to require 
in each case that the term in the error estimate which is dominant in the asymptotic (L — )■ oo) limit is larger than 
the other error terms. This second reason also only requires that L be at least a constant times R. The third reason, 
which in fact gives the most stringent requirement on L and leads to Eq. ([2]), is that our choice of r in corollary ([s]) 
requires an upper bound on r, which depends on L as given in Eq. ( 149| and so requires a lower bound on 07, and 
hence a lower bound on Gr^j^j^L). For simplicity, we summarize all these requirements as "for sufficientUy large". 

While the final estimate ([ij relies on many estimates along the way, we can give a very brief explanation of where 
the scaling with L^^^ arises: in the section on Lieb- Robinson bounds for interactions S^, we are able to derive 
meaningful bounds for a scaling as L^^^ (up to log factors), and then this scaling of a determines the resulting error 
estimates, since the dominant terms elsewhere scale exponentially in the square of a. After giving the proof, we give 
various extensions. First, we give a corollary, using the Hall conductance as an obstruction to finding certain paths 
in parameter space connecting certain Hamiltonians. We then discuss extensions to the fractional quantum Hall case, 
in which we need an additional assumption discussed below which describes the presence of topological order; in the 
fractional case, there are multiple ground states, in which case the phases need not exactly cancel for evolution around 
the large loop. The assumption of topological order deals with this issue, as discussed in the section on the Fractional 
Hall Effect. 

INTRODUCING THE TWISTS 

The evolution we are interested in will be generated by Hamiltonians of the form 

H{ex,(j)x,ey, (j)y) = ^y^- (3) 

ZCT 



(0,0) 



x=L 



FIG. 1: Lines illustrating how the twists are defined on the torus. The twists (jjx, 4>y affect interactions close to the vertical and 
horizontal dashed lines, respectively, while the twists Ox,Oy afTect interactions close to the vertical and horizontal solid lines. 

In order to specify the interaction terms in the Hamiltonian with twists we introduce, for any operator A, the 
following periodic operator twists: 

Qx= 1^ (4) 

l<x{s)<L/2 

Qy= J2 9- (5) 

l<y(s)<L/2 

Then, the flux-twisted interaction terms ^{Z; 6^, (j)^, 9y, (py) are defined according to the following prescription: 
X-1. If 3.se Z : \x{s)-l\ < R, then $(Z; 0„ 0„ 0^, = (0„ $(Z; 0, 0, 0,,)). 
X-2. If 3s e Z : \x{s) - {L/2 + \)\ < R, then $(Z; 0„ 0„ 0^) = Rx{~^x,HZ;O,O,0y,<j)y)). 
X-3. Otherwise, ^{Z;e^,(j)^,ey,(l)y) = <i>{Z;0,0,9y,(t)y). 
Continuing, the terms ^{Z;0,0,9y,(py) are defined as follows: 
Y-1. liBseZ : \y{s) ~ l\ < R, then <i>{Z;0,0,ey,^y) = i?y(6ly, $(Z)). 
Y-2. liBseZ: \y{s)-{L/2 + l)\ < R, then $(Z; 0, 0, = i?y $(^)). 
Y-3. Otherwise, <^>{Z;0,0,9y,(l)y) = $(Z). 

Note that since the interactions ^{Z) commute with the total charge Q and have finite range R, the twists only affect 
2 horizontal and 2 vertical strips, each of width 2R, centered on the lines a; = l,a; = L/2 + l and y — l,y = L/2 + 1 
(the boundaries of the sets on which Qx and Qy are supported), as shown in Fig. [T] Moreover, for each ^{Z) at 
most two "perpendicular" twists act non-trivially to produce the corresponding twisted interaction. In particular, 
there are 4 sets of interactions that see both twists simultaneously, namely those that satisfy conditions X-1 A Y-1, 
X-1 A Y-2, X-2 A Y-1 or X-2 A Y-2. For example, interactions ^{Z) that satisfy conditions X-1 A Y-1 transform into 
RyiOy, Rx{Sx, ^{Z))) (the order in which we apply the twists is irrelevant, since [Qx,Qy] — 0). 

We will need bounds later on the partial derivatives of the Hamiltonian with respect to the twists. There are 
L sites in T (the lines x = 1 and y = 1, respectively) on which the individual terms in d0H{6,O,6y,4>y)\g^g^ and 
dgH{9x, (t>x, 0)\g^g^ may be grouped to act. Each set of interactions acting on a site s has norm bounded as follows: 

J2 \\d9HZ;e,o,ey,^y)\g^gj = J2 \\[Qznx,HZ;ex,o,ey,cj)y)]\\ < q^,^\znx\ J2 Mz)\\ < Qm..J, (6) 

Z3s Z3s Z3s 

where Qznx is the charge in Z that is also contained in Qx, and we defined: 

Qmax = Qmax 



Rx{e,,A) = e^e==Qx^e"^«^Qx ^ i?x(^. + 27r, A), 
RyiOy^A) = e^vQ^ Ae-''^-"'^^ = i?y(6'y + 27r,A), 
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We also used that A > =^ \\[A,B]\\ = \\[A - {\\A\\/2)1, B]\\ < 2p - (|1 A|l/2)1|1 = \\A\\\\B\\ and \Z n X\ < 
\Z\ < R^/A. Similarly, Y.Z3s \\de<^{Z-e^,(t>^,e,Q)\e=eJ < Qmax J- The previous bounds imply: 

\\dgH{e,Q,ey,^y)\e^eA\ < Qm..JL (8) 
\\dgH{0,,<j>,,0,O)\e=eJ\ < Qm..JL. (9) 

Finally, remembering that commutes with the total charge Q, one may verify using the definitions of twisted 

interactions, that: 

Rx{-e, H{B.,,4>.,,By, ^y)) - H{e, -e,<t>^ + e, ey,^y) (lo) 

H{0.,,<j>.,,0y, Cf>y)) = H{9,, 0„ 0y -9,^y+ 9) (U) 

which implies that the twist/anti-twist Hamiltonians (i.e. (p^ — —9x and (j)y ~ —dy) are unitarily equivalent to the 
original gapped Hamiltonian. In particular, we have: 

Rxi&x, RviSy, Ho)) — H{9x, —9^, 9y, —9y). (12) 



THE QUASI-ADIABATIC EVOLUTION 

An important ingredient of cyclic adiabatic evolution is the persistence of the spectral gap in the Hamiltonians 
driving the evolution. The quasi-adiabatic evolution defined below will approximate the adiabatic evolution so long 
as the gap is sufficiently large. Often times, though, we simply don't know how the spectral gap behaves during the 
evolution, but using virtual fluxes we will be able to show that a state quasi-adiabatically continued around certain 
cycles has energy that is very close to the ground state energy of the original gapped Hamiltonian. This observation 
is enough to allow us to deduce that the evolved state is close to the initial non-degenerate ground state, up to a 
phase factor. 



Approximating the adiabatic evolution 

In order to approximate the adiabatic evolution of the ground state of a family of gapped Hamiltonians, we first 
need to write down a differential equation describing the desired evolution. To do this, we note that if for some s > 0, 
we have that {H{9)}g^[o^s] is a differentiable family of gapped Hamiltonians with a differentiable family of unique (up 
to a phase) ground states |^o(^)) and ground state energies Eo{9), then differentiating {H{9) — Eq{9)) \^q{9)) = 
yields: 

(1 - P,{9))de \^,{9)) = doH{9) \^,{9)) , (13) 

where Po{9) = |^'o(6')) (^'o(0) |. Since the phase of \'i'o{9)) is arbitrary, we may choose 

\^'oi9)) = e~IoC^o{t)A'S'o(.t))dt \^,^(^g)) . (14) 
which is also a ground state of H(9) and satisfies the parallel transport condition: 

{%i9),d9%i9))^0. (15) 



Combining the above condition with ( 13 1 we have 



de K{9)) = ^ hIq^^HIo) do H{9) \^',{9)) . (16) 

From here on, we will assume that whenever we adiabatically transport a unique ground state along a path in parameter 
space, we do this to satisfy the parallel transport condition. Moreover, if H(9) has a spectral gap A(9) > 0, then we 
have A{9) ■ (1 - Po{9)) < H{9) - £'0(6') which combined with ^ imphes the bound: 



\\de me)) II < (17) 
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Moreover, we get after expanding dg[{H{9) — Eq{6)) |^'o(6'))] = and rearranging terms: 

dl \%{e)) + {{de%{9),de%{9))) \%{e)) = [2deH{9) {de \%{9))) + {dlH{9)) \%{9))] 



m 

i\\deH{9)r WdlHm 
- A2(0) + A(0) ^^^^ 



where we used condition ^ to get {^'f^{e),d'^g%{9)) = - {dB%{9) , dg%{9)) from de {^'f^{e),dg%{9)) = 0, for the 
second hne and used (17 1 to get the final estimate. 

It will be useful to have an estimate on the gap A(6') as 9 changes from 0. 

Lemma 1. IJ for some s>0, we have that {^(^)}ee[o.s] ^ dijjerentiable path of Hamiltonians each with spectral 
gap A(9) and ground state |^'o(0)), then for 9 £ [0, s]: 

A(6') > A(0) - 2|6i| • sup \\drH{r)\\ (19) 
re[o,e] 

Proof. Let Eq(9) be the ground state energy of H{9). The spectrum of H{Q) is contained in {_Eo(0)} U [-Eo(O) + 
A(0),oo). Therefore, a well-known result in functional analysis implies that the spectrum of H(9) is contained in 
[£;o(0) - \\H{e) - H{Q)\\,Eo{Q) + \\H{9) ~ iJ(0)||] U [Eo{{)) + A(0) - \\H{9) - i/(0)||, oo). We have \\H{9) - H{Q)\\ < 
\\drHir)\\dr < \9\ ■ sup^^^^S ||a,i/(r)||. Thus, A{9) > A(0) - 2\9\ ■ sup^^fo^,, \\drH{r)\\. □ 

Therefore, the gap A{9) remains non- vanishing as long as \9\ ■ sup^^jo e] \\^v-^i''l)\\ < ^(0)/2- 



The super-operator Sa{H,A) 



Now, combining ideas from [T] and [5], we introduce the following super-operator that will allow us to approximate 
the adiabatic evolution given by (16 1: 



Sa{H,A) 



Sait) 



T^{A)du 



dt. 



(20) 



with T^{A) = e'^Hj^^-zuH 



and 



Sa{t) 



V2Tr 



e 2° 



a 



(21) 



The free parameter a can be fine tuned later to give the best approximation possible. Note that the following useful 
bound holds: 



\\S^{H,A)\\<2j^ s^{t)i^JJ\T:^{A)\\dujdt^2\\A\\ s^{t)tdt ^ -^\\A\\ 



2a 
V2^ 



\Al (22) 



where we made the substitution s = t^ l{2c?) in the previous to last equality. Note also that for Hermitian H and A, 
the operator 5Q(-ff, A) is also Hermitian (which will be essential in error estimates later). 
We define now the following state as an approximation to the true ground state of H{9): 

de\<^c.{0))=iSc.{H{9),deH{9))\<i>^{9)), |$,(0)) - |*o(0)) , (23) 

where |^'o(0)) is the ground state of H{0), assuming that such a state is unique. Defining the unitary Ua{9) by 

dgUa,i9) = i Sa{H{9),deH{9)) U^{9), [/„(0) = 1, (24) 
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we may also write: 

\<^c.i0)) = U^iO) |*o(0)) . (25) 

We say that the state \^a{S)) is the quasi- adiabatic evolution of |^'o(0)) along the simple path 0^9. 

The following lemma gives precise bounds on the deviation of the above approximation from the true ground state 
of H{0): 

Lemma 2 (Adiabatic approximation). Assume that for some s > 0, {^^(^)}ee[o,s] ^•^ differentiable path of gapped 
Hamiltonians, each with spectral gap A(0) > 0. For \<^a{0)) defined in and |^'o(^)) the ground state of H{9) 
satisfying (15), the following bound holds for 9 € [0, s]; 

m^{9))~\Mm<2\9\ sup \\de'H{9')\e,=,\\ . (26) 

»)e[o,e] ^(f) 

Moreover, we have the bound: 

\\tS^{H{0),deH{9)g=o)\^om - del^om \e=o\\ < \\deHme=o\\ ^^^^ ■ (27) 
Proof. The above estimates follow from considering the norm of the vector 

M9)) ^ \^^{9)) - , 

which satisfies |(5q(0)) =0 and 

de\6^{9)) ^ iSa^{9)\Sc,{9))+t{Sc,{9)-T{e))\^o{0)), (28) 

where 

S^{9)^So.{H{e),deH{9)), ^T{9) = -^^^y^^^^ deH{9). (29) 



and the definition of iT{9) follows from (16). 

We note that for {|^'r!(6'))}ri>o eigenbasis of H{9), we have: 



(vl/„(0),z5„(0)vl/o(0)) = i^i J Sait)[J^ e™(^"W-^«W)dujdtj (vl/„(0),a,77(0)vl/o(0)) 

. (/>(.) (.---™ - ..) ^■z;"Zr " ^ 

{^o{0),iSa{9)'i'o{9)} = (i [ s^{t)tdt){^o{0),dgH{9)^o{9))^0, (31) 



where the last equality follows from Sa{t) being an even function. Now, setting 

\^'{9))^ziS^{9)~T{9))\Me)), (32) 

we have that : 

{^nie),nO)) = (/I W e'-^""'^^-""^^»^t) , n>l, (33) 

(*o(e),*'(e)) = 0, (34) 

where En{9) is the eigenvalue of |^'„(6')). Recalling that E„{9) — Eo{9) > A{9), for n > 1, the above equations and 
the Fourier transform of the Gaussian Sait) readily imply: 

\\\^'m\\<\\do'H{9%,^e\\^^, (35) 
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where we used the identity ||^ l^*) |p = X]n>o I (*n(^)' P the norm inequality \\A j^f) \\ < \\A\\ for unit vectors 
l^f). Finally, setting 



and recalling the definition of Ua{0) in (24) we have: 



A«(0) = C/„W ( Ulit) {W{t)){5^{t)\ + |<5„(t))(vl/'(i)|) dt\ ul{e), 



(36) 



which one may verify by differentiating with respect to 9 and then using (28 1. Now, taking the trace on both sides, 
we get: 



\UeW = / {{5o.{t),-^'{t)) + {^'{t),6^{t)))dt 



< 2\9\ sup \\5^{t)\\ sup \\^'{t)\\ 
te[o,e] te[o,e] 

Choosing te e[Q,e]:\\5c,{te)\\ — supjgjQ g] the above bound implies: 



¥o.{teW <2\te\\\5M\\ sup 

t6[0,t8] 



1^.(^^)11 < \\5M\\<2\tg\ sup 

te[0,te] 



(37) 
(38) 

(39) 



Combining the above bound with (35 1, we immediately get (26 1. Now, note that (28 1 combined with (35 1 imply (27 1 
as follows: 



\\do\6^{9))\e=,\\^\\ Wm\\<\\deH{e)\e=,\ 



A(0) ■ 



where we used (22 1 in the estimate above. This completes the proof. 



(40) 

□ 



Quasi-adiabatic unitaries and loop operators 



Since we will be looking at quasi-adiabatic evolutions corresponding to Hamiltonians H{9x,9y) with domain em- 
bedded in (actually, {9^, 9y) G 27r x 27r), we introduce here the unitaries describing the quasi-adiabatic evolution of 
l^o) — 1^0(07 0)) around a closed loop in flux-space. We begin by defining below the generators of the quasi-adiabatic 
dynamics: 



'DxidxjSy) — Sa{H{9x,9y), dgH{9,9y)\e=e^) 
T>Y{9x,9y) = Sa{H{9x,9y), deH{9x,9)\g^0^). 



(41) 
(42) 



Note that for notational purposes we suppress the dependence of these operators on a. We continue with the unitaries 
corresponding to evolution through changes in the 9^ and 9y fluxes, respectively: 



X-) ^y-) 

r) = i'Dx{9x + r,9y)Ux{9x,0y,r), Ux{9x,dy,0) = 1 
drUY{9x,0y,r) = iVY{9,,9y + r)UY{9x,0y,r), C/y (0„ 0) = 1. 

It will be useful to note here the following composition rule for Ux{9x, 0y, r), Uy{9x, (^y,f)'- 

Ux{9x + r,9y,s)Ux{9 r) = Ux{9x,&y,r + s), 
UY{9x,Oy + r,s)UY{9x,0y,r) = UY{9x,9y,r + s), 



(43) 
(44) 



(45) 
(46) 



which is easily verified upon differentiating both sides with respect to s and is equivalent to evolving for "time" r + s 
in the X and Y directions, respectively, by evolving first for "time" r and then for "time" s. 
Using the above unitaries, we construct the following useful evolution operators: 



V[{9^,9y) ^ {(j)^,(j)y)] = Ux{0x,4>y,<t>X-9x)UY{9x,0y,4)y-9y) 

w[{9^,9y)-^ {<l)^,<l)y)] = UY{<t>x.0yAy-ey)Ux{ex:ey,<t>^,-9^) 

Vo{9^,9y,r) = V^[{9^,9y)^{9, + r,9y + r)]W[{9,,9y)^{9, + r,9y + r)] 



(47) 
(48) 
(49) 
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Note that V[{9x,9y), {(j>x,(t)y)] can be thought of as evolving a state along the path Ty ■ {Ox,Oy) — > {d^, 4'y) — > {(j>x, <j>y) 
in parameter space, while Oy)^ {(j)x, 4'y)] would evolve a state along the path Tw : {9x,0y) {(px, 0y) {<t>x, 4>y)- 

Finally, VQ{9x,0y,r) is equivalent to evolving counter-clockwise around a square of side r and origin {9x,0y). 
We introduce now the following important family of states: 

l*o(0> = ^o(0, 0, r) \^o{0, 0)) = 4(0, 0, r) uj,{0, r, r) Uyir, 0, r) f/x(0, 0, r) \^o{0, 0)) (50) 

These states describe the quasi-adiabatic evolution of the initial ground state |^'o(0,0)) around a square of size r, 
starting at the origin in flux-space. The estimate (26 1 directly implies the following phase estimate for an adiabatic 
evolution around a closed path 



A(r) = (0,0) ^ (r,0) ^ (r,r) (0,r) ^ (0,0), 
where the arrows represent straight lines between two points in flux-space: 



(51) 



Corollary 1. Assume that the two-parameter path of differentiable Hamiltonians {H(9x,0y)}(^g^ g )eA(r) maintains a 
uniform lower bound A > on the spectral gap. Then, for |^'o('')) defined in {50) and |^'o(0a;, f?j,)) the ground state 
of H{9x,9y), the following phase estimate holds: 



where 



(vl/o(0,0),vl/o(r))-e*^('-) <8rDH{r) 



dOx ln,{{^o{0x,r),de^-^^{9x.r)) ~ {■^o{9x.Q),de^-^o{9xM} 
d9y Im {(*o(r, 9y),de^-^o{r, 9y)) - («'o(0, flj,), ae„*o(0, 9y))] 
d9x f d9y lm{{dg^^o{0x,9y),de^^oidx,0y))} , 



(52) 



= 2 



DH{r) = sup {\\deH{9,9y)\e=eJ,\\deH{9x,9)\e^eJ} 



(53) 

(54) 
(55) 



and Im{-} stands for the imaginary part of a complex number. 

Proof. Deflne the following four vectors, which correspond to deviations of the quasi-adiabatic evolution at the end 



(56) 
(57) 
(58) 
(59) 



of each leg in A(r): 










|<5c.(r,0)) = C/x(0,0,r)|vI/o(0,0))-e^'^^('- 


l^o(r,0)). 


101 (r) = - 


- / (*o(6 
Jo 


x,0),de^^o{9x,0))d9x 


|<S„(r,r)) = [/y(r,0,r)|*o(r,0)) -e^^^M 


l*o(r-,r)) , 


«02(r) = - 


- r(*o(r- 


,9y),dey-^o{r,9y))d9y 


|5,(0,r)) = C/i(0,r,r)|vI/o(r,r))-e^*^'('-) 


l*o(0,r)). 




f («'o(ex, 




r),de^^ai9x,r))d9x 


|^„(0,0)) = C/|(0,0,r)|vI/o(0,r))-e»'^^M 


l«'o(0,0)), 


i4>i{r) = 


/ (*o(o,^ 

'0 


g,ae„*o(O,0,))d0j,. 


Then, we have the following decomposition: 

\^Q{r)) = e^'t'^^'^)+^'t>^i'^)+i'l'3ir)+i4>i 


M|vI/o(0,0)) 


g«0l ('•)+» 




|'5a(0,0)) 



+ e»*iM+»*2('-)[/|(0, 0, r) |(5„(0, r)) + e**i('')[/1;(0, 0, r) [/^(O, r, r) |(5„(r, r)) 
-I- [7^.(0, 0,r);7l(0,r,r)f/y(r,0,r) |J„(r,0)), 

which implies: 

(vl/o(0,0),vl/o(r)) - e'Mr)+^Mr)+^Mr)+^Mr) < ||<5„(o,0)|| + ||(5„(0,r)|| + P„(r,r)|| + ||5„(r,0)|| 



(60) 



(61) 
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Now, from Stokes' Theorem we have: 

^(|>ir)=^Mr)+^Mr) + iMr)+^Mr)=2^ I dO, i dOy lui{{de^-^o{ex,Oy),dg^-^o{0.,ey))] , (62) 



so it remains to show that the vectors 5a have small norm. But, this follows immediately from (|26| and the fact that 
the phases {(j)i{r)}i<i<i are chosen so that the respective ground states satisfy the parallel transport condition (151. 
In particular, applying (26l successively to the gapped Hamiltonians H{r,0), H{r,r),H{0,r) and H{0,0), we have 
the bounds: 



ll'5a(r,0)|| < 2r sup \\dgHi9,0)\e^^\\ 

))G[0,r] 

|i<5„(r,r)|l < 2r sup \\deHir,0)\g^J 

rie[0,r] 

ll'5a(r,0)|l < 2r sup \\dgH{9,r)\g^„\\ 
||(5„(0,0)|| < 2r sup \\dgH{O,0)\g^J 



A 



A 



which completes the proof. 



(63) 
(64) 
(65) 

(66) 

□ 



Note that the quantity Im {(^dg^'i'o{0x,0y),dg^'i'o{9x,9y))} is gauge-invariant; in particular, it remains constant 
under phase changes \%{ex,9y)) = e'-'^'^-^f) \^o{ex,ey)), with f^e^^ey) any real, differentiable function of Ox and 9y. 

We now take a closer look at the phase (/)(r) accumulated during the adiabatic evolution of |^'o(0,0)) along the 
closed path A(r) defined in ([sT]). 



Corollary 2. Let (j}{r) be the phase defined in \53 54-). Then, the following hound holds for r > 0: 
|0(r)/r2 - 2 Im { {de^^o{0, 0,)^=o, de^^o{ex, 0)e.=o)}\ 



< 



I 1 r 

sup \ 7.\\[ds9{s,s)]s=s3\\ + Y\ ■ sup \\dl[g{s,S4) + g{si,s) 

S3e[0,r] \ ^ S4e[0,r] 



(67) 



with g{s, s') = 2Im{(9,-.«'o(s, s")s"=s', 9s"«'o(s", s')s"=s)}. 



Proof. The above bound follows from the Taylor expansion of ^(j') up to order 3. We have from (54) that 0(0) = 

dr(t>{r)\r=Q — and 9^0(r)|r=o — 4Im |(9ej^5'o(0, 6*^), 90^^'o(0a;, 0))^ ^^j. This follows from the Taylor expansion 

of the general integral operator K{r) — dOx Jq dOy g{6x, Oy), where g{Ox^ Oy) is doubly differentiable for Ox, Oy G [0, r]. 
In particular, we have 



[dlK{s)\ = 3[5,5(s, s%^s, + r dsi [d!{g{s, s^) + g{s,, s))l 

Jo 



(68) 



which follows from higher order partials of {dsK{s))^^^^ — J^^ ds[g{s,si) + g{si,s)]. Then, expanding around ?' = 
up to third order, we get: 



K{r) = .g(0, 0) + / ds_ 
Jq 

The above expansion implies 

\K{r)/r^ ~ g{0,0)\ < 



dso 



dsrt 



3{dsg{s,s))s=s3 + / ds4,[d'^ig{s,s4,) + gis4,s))]. 



sup 

S3e[0,r] 



\[9sgis,s)]. 



4! 



sup \\d'^[9is,Si) + g{s4,s)]s=sj 

S46[0,r] 



(69) 



(70) 



In our case, g{s,s') = 2Im {(9s"^'o(s, s")s"=s'7 f^s"^o(s", s')s"=s>}- 



□ 
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LOCALIZING THE OPERATOR Sc{H,A) 



We turn our focus in showing that the operator Sa{H, Az), with H = J^xcT^i-^) G Az, is effectively 

localized around Z; the only requirement for the interactions ^{X) of the Haniiltonian H is that they have finite 
range R and finite strength J. For Af > R, let us write: 



Hm{Z) = J2 *(^)' ^it^ ^(^^) = {X CT: d{X, Z) < M ~ R}, 

X<ZZ(M) 



(71) 



for the interactions restricted to within a radius M of the set Z . Then, for any operator Az with support on Z d T, 
we define: 



AI{H,Az) = t^^'^^HAz)~t^''-^^^HAz), ioYk>R, Al{H,Az)=Az. 



(72) 



Note that the action of A^{H, Az) depends on the support Z of its second argument, but we suppress this dependence 
for notational purposes. It follows that for R < M < L we have: 



M 



r^^''^^\Az) = Y.Al{H,Az). 



k=R 



From the above decomposition, we have for Az G Az 

M 



/oo i^t 



It will also be useful to define the following operator here for A = J^zct ^z- 

Sj''HH,A)^Y.^-^"^iH,Az). 



dt (73) 



(74) 



ZCT 



Now, for any Az £ Az with diam(Z) < R and Am{Z) — Hm+i{Z) — Hm{Z), the following Lieb- Robinson bound 
holds: 

rf^^+i(^'(^z),AM(^)l < 2\Z\ ■ ||AAf(Z)||||Az||e''"l"l-'^"^(^^'^*^(^» < Cfl„7|| Az|| Me^l^l-^^/^ (75) 

where C^^j — 2R?J and v — 132e (1 + i?)'^ J, in the last bound. The first inequality, known as a Lieb- Robinson bound, 
and the constants involved are obtained in [7|. The second bound follows from the observations below: 



1. The distance between the terms Az and Am{Z) is calculated using the definition of Hm{Z) in (71): 

d{Az,AM{Z)) = d{Az, Hm+i{Z) - Hm{Z)) = M - R. 

2. For Hamiltonians embedded in the torus T, we have the bound: 

\Z\ ■ \\Hm+i{,Z) - Hm{Z)\\ < {R/2f ■ {AJM) < R^JM. 



3. The Lieb- Robinson velocity v is given by = 2Co(l)||«'||?, for a = 1/i?, with Co(l) defined in ([93]) and 

defined in (94), with /(si,S2) = J2zbsi s2 11*^(^)11 - bound ||<i>||° by choosing d(si,S2) = R, since each 

interaction ^(Z) has range at most R. 



Moreover, using the above Lieb-Robinson bound, we have the estimate: 



Al,^,iH,Az)\\ - rd,[r":'^^^{r^^^^^(^\Az)))^ 



dt 



< 



\u\ 



Hm + i{Z) 



{Az),Hm+i{Z)-Hm{Z) 



dt 



< CR,j\\Az\\Me-^'^'^ r e"l*ldi 

< -i-||Az||Me-*^/«+"H. 
ZR 
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The equality in the first hne of the above equation follows because the integral over t is equal to 
^Hm(z) ^^Hm+i{Z)^^^^^^ _ ^^.^^ ^^i^-j^rily equivalent to Tu"^'^^\Az) - t""^^\Az). The inequaUty in the 

second line follows from 9, 



(r5r*^^ {t^"^'^^\Az))) = it":'^^^ ([Hm+i{Z) - Hm{Z),t^"^\Az)\) and using a 

d: 

<nAz\\, (76) 



triangle inequality. Of course, for all M > R we also have the naive upper bound: 



\Al,+,{H,Az)\\ < 



{Az) 



{Az) 



which is useful for larger times u. 

We now estimate the error we make by using Sa '^'^\H, Az), defined in (73), to effectively localize the operator 
Sa{H, Az) to within a radius M of the support Az- 

Lemma 3 (Localization of Sa{H, A)). For any operator Az G Az, with diam(Z) < R, and H a Hamiltonian with 
finite range R and finite strength J interactions, the following bound holds: 



\\S^{H,Az)~SJ''\h,Az)\\< \\Sc.''+\H,Az)\\ < ^\\Az\\g^{M/R), 

N>M 

where for a = 2av > \p2/n: 

'$,R\f^ae^i^ for 1 < X < a'^ 
for X > 



9a{x) 



32R 



xe 2 



Proof. We begin by bounding the norm of Sa+'{H,Az), for N > M, as defined in ([73|, smce 

L 

Sc,{H,Az)-SJ^'\h,Az)= Sc.^^\H,Az). 



N>M 



Using the estimates on || A]^^-^(iJ, A^)!! derived above, and noting that Sa{t) < — 75=, £ M, we have 



\\S^''+\H,Az)\\ < 2 s^{t)(^j\\Al^,{H,A, 



■)\\du dt 



1 



< 2{2\\Az\\) / So.{t)tdt + 2—\\Az\\Ne-''^'' 



1 







a\/2n 



4a . , 
< -^\Az 



V2tt 
4a 



(T„/a)V2 



e^'dt+ ^= Az -j.—e 
V2tt 2Rav 



2-n ^PtK "2i?a2w2 



N/R 
N/R+vTn 



e^'dt 







4a 



< 



\/27r 
4a 
2^ 



\Az\ 
\Az\ 



1 



1 



N 



AR oP-n" 



e{N/R)i 



1 + cr- 



R 



,-£(JV/_R)f 



where we set 



a — 2av, e{x) = 1 



(77) 



(78) 

(79) 
(80) 
(81) 
(82) 
(83) 
(84) 

(85) 



after choosing Tn so that {T-^ j ct)^ (2 + wT/v — N jR = 0. We note here that t{x) is a positive, increasing function of 
X > 0, with the following properties: 



x> 



<i(x) > 1/2. 



(86) 



Moreover, 
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Using the above observations, (84 1 implies for M > a^R > 2ttR: 



11^ ;v+i/„ < Ml 4a , v-^ 2 iV « 16a , ,, 4i? / A/' 
J2 \\Sc.''+\H,Az)\\ < \\Az\\ -^u^"™ ^ -J^r ll^^ll ^ I ^ I e"™' 
N>M ^27r K ^2tt 



a2 V R 



where, noting that Ne is decreasing for N > 2R, we used a integral approximation for the following sum: 

w>A/ \ N>o N>o y 1-e ™ V ^ 7.= ^ 

<4{MR + R^)e-^ <8i?2 (^^)e-^ft. 

To get the last inequality we used the bound 1 — > x/2 for < a; < ln2, noting that R> 1 => l/(2i?) < In 2. 
For R< M < a^R: 



N>M 



< 



< 



V2Tr 

4a 

V2^ 
16a 
V2^ 



Ra^ 



\A, 



^(2 + V2^aR^ e^^^ + 16i?e" 



R\\Az\\ V27rcr e^^^ 



(89) 



where we assumed <j > \/2tt to get the last inequality. Combining all of the previous observations with the bounds 



(881 and (89), we get the desired behavior of ga{x), as claimed in the statement of the Lemma. 



□ 



Lieb- Robinson bounds for interactions Sa''{H,A) 



We are interested in Li eb- Robinson bounds for systems evolving under Hamiltonians with interaction terms 
Sa'^iH, Az) defined in (73 1, where k > R and Az € Az- We note that each term has support on Z{k), the set 



of radius k centered on Z, defined in (71 1. It will be useful to define the set of supports for all non-zero interactions: 

e = {Z C T : \\Az\\ ^ 0} (90) 

For sites si, S2 G T, setting d{Z; si, §2) — maxi=i.2 d{Z, Si) and S{M; si, S2) = {Z e £ : d{Z; si, S2) — M}, we have: 

5'(M;si,S2) 7^ =^ M> d(si,S2)/2, Z{k) 3 si,S2 =^ k> d{Z-si,S2), (91) 

where the first implication follows from 2M > d{si, Z) + d(s2, Z) and the second one follows from its contrapositive 
statement. Finally, we note that for all si, S2, the sets S{M] si, S2) and S{M'; si, S2) are disjoint for M ^ M' and 



U 5(Af;si,S2) =1 

Af>d(si,S2)/2 



(92) 



For a quantum system embedded in T, the Lieb-Robinson velocity is bounded above by 2Co(e)||$||° [7], where for 
e < 1: 



1 + d{si,S2) 



s'^FeT^^\{l + d{si,z)){l + d{z,.2)) 



2+e 



<34+^ 
e 



and for a > 0: 



||ci>||:= sup (l-fd(si,S2))'+^e'^''(^-^^)/(si,52), 

Sl,S2€T 



(93) 



(94) 
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with 



k'>d( Z-.si ,32) 



M>d{si,S2)/2 \zeS{M -.31,32) 



[\\S^iH,Az)l\\Sa{H,Az)-Sj''\H,Az)\\}^ 



(95) 
(96) 



where we used (91 1 and (92 1 to get the second hne 



The bounds for Co(e) are derived as follows: 



2+e 



sr^ f l + rf(si,S2) y ^2 

Ui + d(si,z))(i + d(z,s2))y - 



Z^Si,32 



? Vl + d(si 



1 + d{si,S2) 



< 2+ E 



Z#Si,S2 



d(si, S2) + d(si, z) • d{z, S2) 

2+e 



2+e 



min{d(si,z),d(z, S2)} 



< 2 + E (^2 m) • (2/m) 
1 



2+e 



m>l 



< 2 + 2 



4+e 



1 



(97) 

(98) 
(99) 

(100) 



where C2 = 4 is the surface area of the unit ball in Z^, the first inequality follows from setting y 
miii{d(si, z),d(s2, z)}, a; = max{c?(si, z), c?(s2, z)}, w — d{si, S2) and noting that: 



{y <x) A {w <2x) =^ {l + w- 2x){y - 2) < 4x =^ (1 + w)y < 2(1 + w + xy) 



-1±^<'- (101) 
1 + w -\- xy y 



and the last bound is derived from X]m>i ™ ^ 1 + ^ dx. 

Now, in order to bound ||tf>||e, we estimate an upper bound for /(si,S2) using (22l and Lemmajsjas follows: 



/(S1,S2) < 



2a 
\f2^ 



J2 mm{l, g^{M/R)} ^ 



M>d{si,32)/2 



2a 

< I sup 

V27r I A/>o 



E 



^^ZeS(A/;si,S2) 

E min{l,5„(M/i?)} 



2a 

< ^=J 
2tt 



ZeS(M;3i,S2) 

max ' -^max ^ ^ 

M>d{si,32)/2 



M>d(3i,32)/2 

min{l,5r„(A//i?)}, 



where we set 



Jmax = sup 

\ SET 



YUz\\] 

Z33 I 



, Smax sup sup |S'(A/;Si,S2)| ) , 

s\,S2ET \M>0 



(102) 

(103) 
(104) 

(105) 



with 15*1 denoting the number of sites in the set S (considered as a union of the sites within its subsets). We turn 
our attention to ha(d) = X]A/>d/2 "which we estimate in three separate cases, using the definition 

of gaix) in Lemma [s] and the fact that 1 < ga{x) for x < ■y/2 ln(8i?cr^)fT. Setting 



we have: 



do = v/21n(8cr2i?) {2aR), 

' (do - d)/2 + iirR for < d < do 
ha{d) < I 327rcr2i?2 ^ 2(2^r)2 foj. do < d < 2a^R 



(106) 
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where we used integral approximations to bound ha{d), as in the proof of Lemma [sj Putting everything together, we 
have: 

ll^lla < ^Jmax-W Hiax {{1 + df+' e'^'' h^{d)} (107) 
VZTT 0<d<L 



Differentiating with respect to d the expression in brackets after choosing the appropriate bound for ha{d), we obtain 
the following uniform upper bound for the given value of the parameter a: 



^rnax^ {(1 + df+^ e'"^ K{d)] < dl+' with a = ^. 



Finally, setting e = (In do) ^, we get: 



2e , 

adl (108) 

V 27r 



Applying the bound in [7j to a unitary U{9) satisfying: 

d'eU{e')e=e' ^iY.So,{H{e),de<f{Z;e))U{9), U{0) = I, (109) 



zei 



with jmax ^ Qmax using and Sjiiax 5: 16i?,^, from the geometry of the set £, we have shown that for a — dg^ 
and e = (hido)~^, with do given in (1061: 



2Co(e)||$||: <4Q„axi?"'rfoyiM^: (HO) 
where we used 2Co(e) < 132/e, for e < 1, to get 

2Co(e) jmax • Smaxa^ < 4(3niax ^"^^0 011 C^O, 

recalhng that v = 132e(l + i?)^J and noting that ^^/lndo < y/2 ln(8cr2i?). Finally, for A e yljc, -B e and 
C T, we get: 

||[f/(0)A[/t(0),i?]|| <2||A||||i3||nM|X|,|y|}e-^+"('?— (111) 



where do is given in ( 106 ) and we define: 

t^(Qmax, do) = 4gmaxi?^^do\/lndo- (112) 



Application to evolutions based on UxiO,Sy,dx 



We end this section with an important lemma relating localized versions of the unitaries Ux{0,dy,dx) and 
Ux{0,0,dx) defined in (43 1 with H{9x,0y) = H{9x,0,0y,0). First, we introduce the unitaries UQ{9x,9y,0) satis- 
fying V(?„0,, e [0,27r]: 

de'Uniex,0y,0')e'=0 - -iSj''/^^\Hiex - e,0y),dg,Hn{0\ey)g,^e^^e)UniOx,0y,0), Un{0.,0y,O) = 1, (113) 
with Hii{0x,0y) = Ezcn^iZ;Sx,O,0y,O) and 

(114) 



n=^seT:\y{s)\<-L-R 



The following Hamiltonian decomposition will be useful: 



H{0x,0y) — HQ{0x,0y) + Hn':{0 )- ^ <f>(Z;0„O,0„O). 



(115) 



Finally, since the unitary Un{0x, 0y, 0) acts trivially outside of the set fly defined in ( 120 1 below, the following is true: 

RYi0y, Uni0x,0, 0)) = Uni0x, 0y, 0) (116) 

The Lieb-Robinson bound we have obtained on the velocity of propagation for interactions of the form Sa^{H, Az), 
has a direct application to estimating the following important upper bounds: 
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Lemma 4 (Twisting L emm a). For UxiO,Oy,Ox), Ux{O,O,0x) defined in (43) with H{0x,0y) = H{9x,0,9y,0) and 
Un{(^x,()y,S) defined in (113), we have for sufficiently large L and a — ( 48 iS L I ^ ' C ^ i^/"^) 



||[/i(0, 0y, Ox) An, f/x(0, Oy^e.,) - Un{,0,,ey,e,) An, uUe.,, By, e,)\\ < |0,|^g„,ax J\\Ano\\ L'e-(j^)''' (117) 

\/2n 



||[/i(o, 0, 0x) An„ Ux{o, 0, e,) ~ Uniex,o, Ox) An, ul{e^,Q, e^)\\ < I^^I^Q^ax 

V 27r 

for all 9x,0y G [0, 2tt] and An, G An, with 

no = {seT: \x{s)\ <L/8-R and \y{s)\ < L/8 - R} . 



J\\Ano\\L^e-(^) , (118) 



(119) 



Moreover, for all 0x,9y G [0, 27r] the evolved operator Un{Oxi(^yi(^x) An, U^{9x,9y,9x) has support strictly within the 
set fix n fly , where: 



nx = {seT: \x{s)\ < L/4}, fty = {s e T : |y(s)| < L/4}. 



(120) 



Proof. Wc only show the first bound, since the second one follows from a similar argument. First, note that 
Un{9x,9y, 9) satisfies: 

de'Un{9x,Sy,9x — 9')gi=e = —dg"Un{9x, 9y, 9")e"=e^~0 = iSa^^^^^''{H{9, 0y) , dg' Hn{9' , 9y)g>=e) Un{9x, Oy, 9^ — 9), 
where we made the substitution 9" — 9^ — 9' to get the first equality and used ( |113[ ) to get the second equality. Set 

Ac/(0x, Qy) = ul,{Q,9y,9,)An,Ux{Q, ) — Un{9x, Oy, 9x) An, ul^{9 X 5 ^y 7 ^x ) 



and 



for 



As{9, 9y) = S^{H{9, 9y),de'Hi9', 9y)e'=e) - S^^'^'^^\h{9, 9y),dg,Hn{9' , 9y)e.=B) 
= Af/'^\9, 9y)+Sa. (^/24) (^(^^ ^ Q^^jj^^ (g,^ 9y)o'^e), 

Af'^^\9,9y) = Sa,{H{9,9y),de'Hn{9',9y)e'=e)~Sj'^/^^HH{9,9y),dg,Hn{9\9y)g,=e) 

+ S^{H{9, 9y),de'Hn^{9', 9y)g,=g) - SJ'^'^^\h{9, 9y),dg,Hn40' , Oy)g'=e). 



Now, we have 



\^uidx,9y)\\ 



de> (ul,{0, 9y,9') Un{9x, 0y, 0x - 9') An, ul{9x,9y,9, - 9') Ux{0, 9y, 9') 



d9 



< \0x\ sup \\[Un{9,,9y,9)An,ul{9,,9y,9),Asi9,9y 



< 2\9,\\\AnJ sup 



AL/24) 



(O.Oy) 



\9x\ sup [Uni9x,9y, 9) An, Ul{9x, 9y, 9),sJ''/^^\h{9, 0),de,Hn^{9' , 0)g,=e)] 
0e[o.e,] 



(121) 



Applying the Lieb- Robinson bound ( 111 I, we have for G [0, 27r]: 



Un{9x,9y, 9) An, ul{9x, 9y, 9),sJ'^/^^Hh{9, Q),deHn^{9, 0)) 



< 



4a 



27r 



(122) 



where do is given in (106), ^(Qmaxi c?o) is defined in (1121 and we made use of (22 1 with the bound ([6]), to get: 



\S^'^'^'^^\H{9,Q),dgHn^{9,m\ < ^Qm..JL. 

V 27r 
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Moreover, |rio| < L'^/GA. Now, setting 



i?2 



4(27r)Q„ 



(123) 



and assuming do < ( 43^! l I ' ™pli'3s the bounds > 27r ^(Qmax, do) and, hence, 5^ - 27r w((3i„ax, do) > 4^ > 



24(io 



(s)^^^ {^^^y^^ ^ (iIr)''''^ foi" t^^c exponent in (122 1. Going back to (122), we get 



Un{o,,ey,e) An,, ul{e^, 9y, 0),sJ'^/^^Hh{9, 0),deHn4e, 0)) 



< 



16V27r 



8i?. ^ 



Finally, using Lemma [s] after noting that > (do/2_R)^ > cr^, we have the bound 



2\\A 



Ho II sup 

9e[o,e^ 



(124) 



(125) 



Noting that ( 124 ) dominates ( 125 1 in (1211, completes the proof of the first bound in the Lemma. To justify the choice 
of the global parameter a in the statement of this Lemma, we recall that we have assumed dg = 2cri?-\/ln(8a^i?) < 



/InL 



and since a = 2au, we choose 



1 



4wi? V481n'^L 



to satisfy the bound on dg while maximizing a. 



(126) 

□ 



COMPUTING THE HALL CONDUCTANCE 



We use the Kubo formula from linear response theory as applied to the setting of a torus pierced by two solenoids 
carrying magnetic fluxes 6^ and 9y in the x and y directions, respective ly|T3]. We denote the Hall conductance at the 
origin in flux-space by a^y, given by the formula: 



a^y = 2 Im 



80. 



09, 



2n- 



(127) 



where 1^*0(0, 0y)) and |^'o(6'a:, 0)) are the ground states of H{0,0,6y,0) and H{9x, 0,0,0), respectively. We note here 
that the flux angles 9x and 9y measure the flux in units of fi/e, the unit of magnetic flux divided by 2tt, so that 
H{9x,9y) is periodic in both angles with period 27r. With this in mind, we may also write for dimensionless 9x,9y: 



To show quantization of a^y, it is sufficient to prove that 



\ _ g2iri CTa;^/(e /h) 



decays as a stretched exponential in L. It follows then, that 3 n e N such that 

c2 



< 



V2e2 

2TTh 



;*o,*o(?-)) 



+ |l-(vl/o,v|/o(2^))| + 



is bounded above by a function that 
(vl/o,vl/o(2^))-(vl/o,vl/o(r))(^)' 



(128) 



where we used 



27r/i „ 
'^xy-^ 27rn 



< V2 



1-e 



i (cr^y 2irn) 



which follows from assuming 1 1 — * I < 1 =4> |1 — e*^p = 2(1 — cos0) < 1 16*1 S [0, 7r/3], up to integer multiples 

of 27r. In that case, cos \9\ is monotonically decreasing from 1 to 1/2 in that range. In particular, we have 



rW r't> r 
= 2(1 -cos|6l|) = 2 / d(/) / cosT]dr]>2 

Jo Jo Jo 



cos77d?7>2 / (j) cos (j) d(j)>\9\'^ cos \9\>\9\'^/2. 



Since |l -e'l^N = |l - e'^l, we get \9\<V2\l-e 
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At this point, we are ready to derive bounds for each of the three terms in (1281. The first one emerges from the 
close relationship between the Hall conductance and the geometric/Berry's phase 4>{r) that appears in Corollary [T] 
In particular, we know from Corollary [l] that for H{9x, Oy) — H{9x,0, Oy, 0) and small enough r, the following is true: 



< 



A 



A = 7 - 2Q„-,^^JLr, 



(129) 



where we used (8|9) to bound Dnir) in (52 1 and the gap estimate (19 1 to compute the uniform lower bound for the 
spectral gap of the family of Hamiltonians {H{6x,0,6y,0)}(^g^^g )gA(r), remembering that if(0, 0,0,0) — Hq has gap 
7 > 0, by assumption. From the definition of A above, we can see exactly how small r needs to be in order for the 
gap to remain open during the quasi-adiabatic evolution. At this point, we note the following useful inequality: 



g4ir irj>{r)/r 



r • A 



A = 7 - 2Qmax JLr. 



To show ( 130 1, we use the general inequality for real b with < 5 < 1 and m € N derived in the Appendix: 

Ib-e'"] < e < 1/2 =^ |6"-e 



< 



The above relation implies that < 6 < 1: 



be" 



pi(e-0) 



< e < 1/2 



■ m _ im{0-<p) 



< 



me. 



(130) 



(131) 



(132) 



which combines with (1291 to give (130). Moreover, combining the simple inequality 

-i{^l-4>2)/2 



gi(01-02)/2 _ g 



2sin|(0i -02)/2| < 101-02 



with (67), (127) and (130), we get: 



Corollary 3. For < r < 



the following bound holds for a numeric constant C > 0: 



< c 



-L 



J 



-L 



(133) 



where we used ( ITpS I and similar bounds for the norms of partials of the form d'^,ds'ds,ds'd'^,d^ acting on 
|^'o(s',s)) and |^o(s,s')); order to bound the partials dsg{s,s) and ds [g{s',s) + g{s,s')] in (67). In particular, 
these bounds are in effect bounds on the norms of the same set of partials acting on H{s' ,0,s,0) and _ff (s, 0, s', 0), 
for which we use arguments similar to those leading to (8][9) to estimate. 

The bound for the second term in ( 128 1 is given as Corollary |4] below. Before we derive the bound, we show certain 
trace inequalities and develop energy estimates that allow us to study the quasi-adiabatic evolution of ^'o(0, 0) under 
a closed path in parameter space, when we can no longer rely on the assumption of a uniform lower bound A > for 
the spectral gap. Note here that we consider each of the paths 

Ai(27r) : (0, 0) ^ (2^, 0), A2(2^) : (2^, 0) ^ (27r, 27r), A3(27r) : (2^, 27r) ^ (0, 27r), A4(27r) : (0, 27r) ^ (0, 0) (134) 

to be closed in parameter space, since Hq = H{2tt, 0,0,0) — iJ(27r, 0, 27r, 0) = H{0,0,2tt,0), which follows directly 
from the Aharonov-Bohm 27r-periodicity of the interactions ^{Z; dx,0, 6y,0). 



Partial trace approximation 

We now prove trace inequalities that describe the partial trace of the quasi-adiabatic evolution of |^'o), both near 
and far from the twists driving the evolution. The estimates after tracing out sites near from the twists driving the 
evolution are based on locality, while the estimates after tracing out sites far from the twists driving the evolution use 
the virtual fiux idea of [T] , by comparing the partial trace to a partial trace of a density matrix in which the angles 4> 
are twisted also. 
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First, we define the following useful states: 



\^x{0)) = Ux{0,0,e)\^o), \^x{0,-0)) = uP{0,0,e)\^o), (135) 
\^y{0)) = UY{o,o,e)\^o), \^Y{0,-9)) = U^^\0,0,e)\^o), (136) 

where the unitaries Ux{0,0,e), U^^\o,0,9), Uy{0,0,9) and ?7^^^(0, 0, 6*) are defined in ( |43p4| ), with dynamics based 
on the one-parameter families of Hamiltonians {H{9' ,0,0, 0)}gi^[Qg], {H{9' , —9' ,0,0)}gi(=[Q^g], {H{0,0,9' ,0)}g>fz[Q^g] 
and {H{0,0,9' ,—9')}gifz[Q^g], respectively. We define the density matrices corresponding to the above pure states by 



Px{9) = l^xmi'i^xm. px{9,-9)^\^x{9,-9)){^SJx{9,-9)\, 
Py{9) = py{9,~9) = \^SJYi9,-9)){^SJY{9,-9)\. 



Now, for the sets Qx, defined in ( |120| , we introduce the following Hamiltonian restrictions: 

Hn^i9) = HZ;9,0,0,0), Hn^J9) ^ ^ ^Z;0, -9,0,0), 

zcnx Z(^nx 

HnA9) = '^{Z;0,0,9,0), Hn^^{9) = ^ ^{Z;0,0,0,-9). 

Using the above restrictions, we have the following useful decompositions: 

H (9, 0,0,0) = Hnxi9)+Hn^^{0), H{9,-9,0,0) ^ HnA9) + Hn^x(9), 
H{0,0,9,0) = HnY{9)+Hn^iO), H{0,0,9,-9) ^ HnY{9) + Hn^(9). 



(137) 
(138) 



(139) 
(140) 



(141) 
(142) 



For any set Z, we define Z to be the complement of Z . We define Q,"^^ to be the set of sites within distance R of 
fix and we define fiy to be the set of sites within distance R of fiy. Then, _ffn^ is supported on Vf^ and similarly 
_ffsi|, is supported on fiy. Thus, iJ is a sum of terms '^{Z) such that each set Z is either supported on fix, or on fi^. 

Lemma 5 (Partial trace approximations). The following partial-trace norm inequalities hold for the states and sets 
defined above, for a numeric constant C > 0: 



max{||Tr^(px(e)-px(0))||i, \\Tr^{pY{9) - PYmWi} < C\9\^ (q^,, 

max{||Tr^(px(e) - Rx{9,Pxm) 111, ||Tr^(py(^) - i?y (0, (0))) ||i} < 

ai L L -2., 



7 / CT 



C\9\ Qina.x — L 
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Proof. We only prove the bounds for px, since the bounds for py follow similarly. Setting H^^\9) — H{9, 0, 0, 0) and 
= H{9, -9, 0,0), we note that 

de'px{9') = i\s^(H''^\9'),dg,HnA9')).Px{9') 



dg.px{9',-9') = I [Sa. [H^^\9'),dg,HnA0')) .Px{9',-9')\ +i[s^ (^H<^^\9'),dg>Hn'=J9')^ ,px{9',~9') 

Furthermore, recalling the definition oi Sc}''^'^\h, A) given in (73 1 and using the fact that dgi Hq,^{9') has support on 
a strip of width 2R centered on the line a; = 1, we have: 

S^i^/^-^){H(^)^e'),dg,HnA0'))=S^^'^'^-''\H(^\9'),dg,HnAe')) 
with support on Moreover, from Lemma [s] with L/4 — R> a'^R, we get: 

16a 



So.{H(^\9'),dg, Hnx {9'))- 5„ ( W4-i^) (^(i) (0^) ^ q^, fjn, {9')) 

So.{H^^\9'),dg,HnA9'))-Sj'^/''-''HH^'H9'),dg,HnA9')) 
So.iH^'\9'),dg,Hn^J9')) - Sj'^/*-''\H^^\9'),dg,HnA9')) 



< 



/2tt 
16a 



,L2 

r2 



'211 



16a ^ __£^^i 
V27r 
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We also note that Sj^/^-^\H^^\0'),de>Hno (6»')) is supported on nx, since dgiHQc (0') has support on a strip of 
width 2R centered on the hne x = L/2 + I. Equipped with the above observations, we turn our attention to proving 
the upper bound on \\Tt^{px{0) — px{0))\\i- 



\\Tr^ipx{0)-pxm\\ 



Tr^[S^{H'^'H9'),dg,Hn, (0')), Px(^')] d0' 



<\e\ sup TT^[Sa.{H<''\9'),de'HnA9')),Pxi9')] ^ \0\ sup sup \Tr{A[S^{H'-'\0'),dg^Hn^{e')), px{0')])\ 



e'e[a,e] 



e'<£[o,e] AeA,,a_ 
11^11=1 



\0\ sup sup \TY{[A,Sa^{H<^^H0'),dg,Hn^{0'))]px{O'))\<\0\ sup sup \\[A,S^{H(-^\0'),dg,Hn^{0'))]\\ 
e'e[o,6»] AEAfjc^ e'e[o,e] Aeyifje^ 

= 1 i|A||=l 

\0\ sup sup \\[A,Sa,{H^'H0'),dg,HnA(^'))~Sj'^^^-''\H^'\0'),de'HnA()'m\ 

e'e[a.e] AeA^o^ 

Ml|=i 



<2\0\ sup \\S^iH^'H0'),dg>HnA(}')) - Sj''^'^''HH^'\0'),de'HnAO'))\\ < \0\^Qn...J 



8'G[0,e] 



2n 



e sj? ' 2 . 



To prove the bound on \\Ti-^{px{0) — Rx{9, Px{0))) \\i, we define the unitary Un^{0) with support on the set 
by the differential equation: 

de' Un, i0') = i 5„ ^^/^-^^ {H^'\0'), dg, Hn, (0') ) Un, {0'), Un, (0) = 1. 

Now, we have: 

AJ0) {px{O)-px{0,-0))UnAO) 



X 



(143) 



\\Tr^{px{0)-px{0,-O))h^ 



Tr; 



Ox 



< 



Tr 



nx 



So.{H('H0'),dg,HnAO'))-Sj'^/^-''HH^'HO'),de'HnAO')),Px{0') 



d0' 



Tr-r 



SAH^^\0'),de,HnAe')) ~ Sj'^'^-''\H'^^\0'),dg>HnAe')).Px{0\-e') 



d0' 



S^{H^'H0'),de,H,rUO')),Px{9',-0') 



d0' 



< 2\9\ sup SAH^'\0'),de,HnAO')) - Sj'^/*-"HH^'\0'),dg,HnAO')) 
e'e[Q,e] 

+ 2\0\ sup SAH^^Hd'),dg,Hnxm ~ Sj'^^''''HH^'H0'),de'HnAO')) 

+ 2\0\ sup SAH^^H0'),dg,Hn^J0')) - Sj''/^-''\H^^\0'),dg,HnA&')) 

96a ^ ^L^ L ,1 

V 27r c 



Finally, noting that px(0, 0) = Px(0) = Pq and that {Rx{0\ px{^i^))}e'e[o.e] is the family of ground states corre- 
sponding to the unitarily equivalent family of Hamiltonians {H(0' , —0', 0, Q)}g'e[Q,6] with spectral gap 7, we bound: 



\\i:r^{px{0,-0) - Rx{e,px{Q,m\\i < \\px{0,-0)-Rx{O,Px{O,Omi < \0\Qn...^Le-""'^" 



(144) 



by applying (26 1 to the quasi-adiabatic evolution of px{0,0). Using the triangle inequality 



\\TY^{px{0)-Rx{O,Pxm)\\i < \\Tr^{px{O)-px{O,-0))\\i + \\Ti-n^{px{O,-0)-Rx{e,Px{O,O)mi 
with the above bounds, completes the proof. □ 
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Energy estimates 



We now prove the energy estimates. 



Proposition 1 (Energy estimate). For the states and |^'y(6')) defined in (135) and (136), respectively, the 

following energy estimate is true, for a numeric constant C > 0; 

\{-^x{e),H{e, 0, 0, 0) -^xm - Eo\ < C\0\ q^,^-l^ 4 + e""'^') 

with Eq the ground state energy of Hq. The same estimate holds for \ {^>y{0), H(0, 0, 9, 0) ^'v'(0)) — Eq\. 

Proof. We will show the bound for \'^x{0)), since the bound for follows from a symmetric argument. Noting 

that Rx{d, px{0)) is the ground state of H{9, —9,0,0), we have: 

EoiO) = TriHoPxm = Tr(i/(0, -9, 0, 0) Rx{9, px(0))) = Tr(i/nx (^) Rxi9, pxiO))) + MHn^„ (O)px(O)) 



where we used the unitary equivalence from (12 1 to get the second equality and also 



TiiH,r(e)Rxie,pxm) = Tr(i7o5,(0)px(0)). 



Now, recalling the decompositions given in (141 ) and noting that px(0) = |^'o)(^o|, we have the following estimates: 

\{^x{0),H{9,0,0,0) ^xm-Eo{0)\ = \Tr {HnAS)iPx{0) - Rx{0, PxiO)))) +Tr {Hn^^mpxiO) - PxiO))) I 

< \Tr{HnAmPx{9)-Rx{9,Pxm)) \ + \Tr {Hn^j0){px{9) - pxiO))) \ 

< \\HnA9)\\ ■ {\\Tv^{pxi9) - Rx{9,pxm)\\i) + \\Hn^jO)\\ ' \\Tv^{px{9) - Px(0))|li 

<JL^ {\\Tr^{px{9) - Rx{e,Pxm)h + \\Tr^{px{9) ~ Px(0))||i) 
and using Lemma [5j completes the proof. □ 

The phase around the big loop. 

In order to derive the next bound, we break the evolution of |^'Q(27r)) around A(27r) into its four individual 
components {Ai(27r)}^^j^ given in (1341 and apply Proposition [l] to each one, which gives: 



Corollary 4. For a numeric constant C > 0, the following bound holds: 



|(vI/o,vI/o(27r))-l|<C|0|Q„ax(^) L3(^^^e-ra+e-""'-), (145) 
Proof. We begin by noting that the 27r periodicity of H{9x, 0, 9y, 0) in each angle, implies: 
(^-0, ^^0(27:)) = (^0, uliO, 0, 2^) [7^(0, 0, 27r) Uy{0, 0, 2^) Ux{0, 0, 27r)^'o^ 



4-^(2^), Po (7^(0, 0, 2Tr) Fq Uy{0, 0, 2^) Pq *x(2^)) + (^'y(2^), Pq Uj^iO, 0, 2^) Qo Uy{0, 0, 2^) Pq *x(2^) 



^^-y (27r), Po [/^(O, 0, 27r) [/y(0, 0, 2^) ^x(27r)^ + (2^), (7^(0, 0, 2^) f/y (0, 0, 27r) ^x{2Tr) 

= I (^'y(2^),^'o) P I {^o,^x{2n)) p + (^'y(2^),«'o) {dx{2n),6Y{2n)) (^-q, ^-jf (2^)) 

+ (^'y(2^),«'o) (^'o,^':v(2^)) ((5;,(27r),(5x(2^)) + («'y(27r),«'o) (^x(2^), t/y(0, 0, 27r) 5x(27r)) 

+ (5y(27r),5^(2^)) («'o,[/y(0,0,2^)^'x(2^)) + (<5y(27r),C/|(0,0,27r),5y(2^)) (^-q, «'x(27r)) 
+ (5y(27r), uliO, 0, 2n) (1 - Pq) C/y(0, 0, 2^) Sx{2Tr)) , 
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FIG. 2: The total evolution around a rectangle of dimension 2 x 3 is decomposed into 6 evolutions that all share the common 
features of evolving from the origin, first with 9y and then with 6^, to reach {6x,9y) in fiux space, then following a tiny counter- 
clockwise loop and finally, reversing the path to go from {6x,dy) to the origin. Any m x n evolution can be decomposed in 
this manner, by completing first the bottom row, as in steps (1) — (3), and then stacking the remaining rows on top, as in 
steps (4) — (6). That is, the unitary corresponding to quasi-adiabatic evolution around the larger loops is exactly equal to the 
product of unitaries corresponding to evolution around the smaller loops. 



where we inserted 1 = Po + Qo between unitaries and set: 

\^x{2n)) = Uxi0,0,2n) I^-q) , \^Yi2n)) = [/y(0,0,2^) I^-q) , 
l^'xi^n)) = C/l(0,0,27r) \^o) , |*Y(27r)) = ^(0,0,2^) \^o) , 

6x{27r) = Qo\^'xi2n)), ,5y(2^) = Qo |^'y(27r)) , S'xi2n) = Qo\^'xi27:)) , 5;-(27r) = Qo |«';.(2^)) . 
The above expression for {'i'o, 'i'Q{2Tr)) combined with the fact that 

I {^Y{27r),^o) P - 1 - ll'5y(27r)||2, | (vI/^(2^), vJ/q) p ^ 1 - ||<5x(27r)f 
and ||<5y(27r)||2 = ||,5;.(27r)||2, \\dx{2TT)\\^ = \\6'x{2tt)\\^, gives the bound: 

I (*o,*o(2^)> - 1| < 2{\\SY{27r)f + ||<5x(2^)||') + 4||<5y(2^)|| ||<5x(27r)|| + ||5y(27r)||2 ||^x(27r)||2 (146) 
Finally, noting that H(2n, 0, 0, 0) = H{0, 0, 2tt, 0) ^ Hq > EqI + 7 Qo, we have 

\\SY{2Tr)f<-\{^'Yi2Tr),Ho^Yi2n))~Eo\, \\6x{27:)\\^ < - \ x{27:), Hq^ xi2n)) - Eo\ (147) 

7 7 

and applying Proposition [T| with = 27r to the above inequalities completes the proof. □ 

DECOMPOSING FLUX-SPACE. 

Now that we have the estimates from Corollary [3] and Corollary |4] it remains to show that the following bound 
holds which is proven below: 

Corollary 5. For a numeric constant C > 0, the following bound holds: 



(vl/o,vl/o(2^))-(vl/o,vl/o(r))(^)^ 



< C {Qn^^^aJ Lf^ (^(QmaxaJL)r5 + <C{Qn,.^aJLf^^e~^, 

(148) 



2tt 



L(Qmaxa^i)'/'el"^'^^J 



(149) 



We proceed with the proof of this bound by turning our focus on a decomposition process used to break the large 
evolution around flux space into evolutions around tiny loops on the (27r) x {2tt) lattice. Fig. [2] describes the process 
used to decompose the evolution around any m x n rectangle into m ■ n evolutions involving the r x r squares which 
form the lattice in flux space. This process is effectively a discrete version of Stokes' Theorem. 
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We define the family of states whose evolution follows the paths that appear in the decomposition process given by 
Fig.|2] 



(150) 



To see why these states are important, note that projecting onto the ground-state after every individual cyclic 
evolution in the decomposition of (^'o, ^'o(27r)) given in Fig. [2] corresponds to the product of the following (27r/r)^ 
terms: 



(*o, *o(0, 27r - r, r)) (*o, ^'c('^, 27r - r, r)) • • • (*o, *o(27r - r, 0, r)) . 



(151) 



The contribution from terms projecting off the ground-state is a bit more involved. It is shown in the Appendix that 
in order to prove the bound (1481 in Corollary [5] it suffices to combine the bound (129) with a bound of the form 
proven below: 



, ye^,ey^%2n] (152) 



Lemma 6. For sufficiently large L, we have: 

I {^>o,^o{e^,ey,r)) - («'o,*o(0> I < C iQm..aJ L)\^ [{Qm.^aJ if + e'"' 
with a given in Lemma^ 

Note that the above bound allows us to translate the coordinates {9x, dy) to the origin (0, 0) in flux space, where we 
can make use of the gap 7 > 0, to get good bounds on the approximation of the evolution of the ground state |\E'o)- 



The Translation Lemma 



We now focus our attention on proving (1521, but before we start, we develop the following important Lemma: 



Lemma 7. Let A^o G -4j2o; with flo defined in ( |119[ ). Then, the following bound holds for a numeric constant C > 0: 

|(^'o,Ao„*o) -(*o,l^^[(0,0) ^ ie,,ey)]RYi0y,Rx{d.,Ano))V[iO,O) i9x,0y)]^o)\ 

J 



Q 



2 2 

max I 

7 



<C{\9,\ + \9y\)\\An„ 
with a given in Lemma^ 

Proof. We recall the density matrices px{9x) and pyiOy), defined in ( |137| - [T38| ). Using the rotation identity (II6I 
for the auxiliary unitary Un{9x,0y,9x) defined in (1131 and recalling the definition of y[(0,0) — ^ {9x,9y)] in (47), we 
get the following equalities: 

(*o, [(0, 0) i9x,9y)] RY{9y,Rxi9x,An„)) V[{0, 0) ^ (0., 0,)]*o) - (*o, An.^o) 
= Tr (pY{9y) C/i(0, 9y, 9x) RY{9y, Rx{9x, An,)) C/x(0, 9y,9,)) - Tr(PoAoo) 



= Tr (^pY{9y) [(7^(0, 9y, 9,) RY{9y, Rx{9,, AnJ) Ux{0, ) — Un{9x,9y, 9x )RY{9y,Rx{9x,An„))Ul{9 

X 1 1 ^x ) 1 

+ Ti{[py{9 

v) RY{9y,PQ)\ Un{9x,9y,9x) RY{9y, Rx{9x, A^o)) U^ii^x^ ^y^^x)} 
+ Tv(^RYi9y,Po)RY {9y, [Uni9x,0,9x)Rx{9x,Ano)uli9x,0,9x)-Uj,i0,0,9x)Rxi9x,An„)Uxi0,0,9x)])) 
+ Tr{[px{9x)-Rx{0x,Po)]Rx{0x,An,)) (153) 



Finally, recalling the sets ilx and fiy defined in (120) and the comment preceding their definition, we have: 



|(*o,^^^[(0,0) ^ {9x,9y)]RY{9y,Rx{9x,Ano))V[{0,0) ^ {9x,9y)]^o) - (*o, A^^^^-o) | 

< uj^{0, 9y, 9x) RY{9y, Rx{9x,Aq„)) Ux{0, 0y,0x) — UQ{9x,9y,0x) RY{0y, Rx{0x, A^g)) ul^{0x,9y,9x) 

Un{0x,0, 9x) Rx{9x, An,) C//,(0x, 0, 9x) - [/^(0, 0, 9x) Rx{9x,An„) [/x(0, 0, 9x) 
+ IIAoJI {\\Tr^ipYi0v) - RY{0y,Po))\\i + \\Ti-n^ipx{0x)-Rx{9x,P,))\\i) 

< WAnoW (Qmax ^ (\9x\4^L' e-i^T'^ + C{\9x\ + \9y\) (^3^ ^ + 



(154) 
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where we used Lemma |4] and Lemma [5j Finally, the condition on a in Lemma [4] implies that the exponential term in 



( 154 1 corresponding to: 

2 2 t 1 ^2/5 / ^ s4/5 

g-a 7 = e"^i^™Mi^i^j >a7L2e"ii5«j ^ (155) 
dominates the other two terms for sufficiently large L. Putting everything together gives the desired bound. □ 

Localizing the loop unitary V,'j{0x,0y,r) 



Now, we show that up to order 4 in r, the operator 1/(3(0, 0,r) defined in (49 1, can be approximated by the sum of 
operators with support on a cross of radius L/8 centered at the origin x = y = Q. Moreover, we show how we may 
approximate V(^{9x^ Oy, r) by Ryi^Oy, Rx{Ox, 1^0(0, 0, r))), up to a small error in depending on L and r. 



Lemma 8. For sufficiently large L, there exists a constant C > 0, such that for all 9^, 0y G [0, 27r] we have the bound: 
\\Vo{e:,,ey,r) - RY{ey,Rx{e:,,Vo{0,0,r)))\\ <c(^ (q„...-l) I^e-^ r + {Q„,,^aJL)\''] . (156) 



vv^27r V 7 / / 
Moreover, there exists an operator W{r) £ Aqo, such that \\W{r) — 1\\ < C {Q^i^y^aJ Vf' r^ and: 

||Fo(0,0,r)-l^(r)|| <C(^^ (^Qmax^i) ^e-5fer + (gmax«Ji)'V5^ . (157) 

Before proving the above Lemma, we introduce the unitaries Ux(M){Qx,Qy,r) and UY{M){Ox,(^y,r) defined by the 
following differential equations: 

(M){(^x,Oy,r), Ux(M){0x,Oy,O) — 1 

UY{M){Ox,dy,0) = 1, 



with iSa'*^-'(-ff, A) defined in (73 1. The following Lemma gives us a bound on the error of approximating the unitaries 



Ux{OxiOy,r) and UY{0x,0y,r) with Ux(M){(^x,Oy,r) and UY(M){Qx,Qy,r), respectively. 
Lemma 9. The following hounds hold for all 6^, Oy G [0, 27r] and r > 0, with M > a'^R: 

\\Ux{ex,Oy,r)-UxiM){Ox,Oy,r)\\ < 64^ (Q^^^^lA^e-ikr, (158) 

||C/y(^^„0J„r)-[/y(M)(^?.,^?„r)|| < 64^ fgmax-i) ^e-rar. (159) 

V ZTT \ 1 J ^ 

Proof. We only prove the bound for Ux{dx, (^y,i"), since the bound for Uy{Ox, dy, r) follows from a symmetric argument. 
First, note that \\Ux{9x,dy,r) - Ux(M){0x,9y,r)\\ = \\Uj^{0x,dy:r)Ux(M){Ox,dy,r) - Moreover, setting 

Aa/(^„ Oy, r) = Sj^'\Hie, + r, 0, 0y,O),dgHie, 0, Oy, O)0^0^+r) - S^{H{e, + r, 0, By, Q),deH{e, 0, ey,Q)e=e^+r) 
we get after differentiating with respect to r: 

dr{y^x{^x.ey,r)UxiM){ex,ey,r))^^^ = i (ul,{9x,9y,s)AM{9x,0y,s)UxiM){Ox,Oy,s)) (160) 

and hence using a triangle inequality we have 

\\Ul,{e.„9y,r)UxiM){0x,Oy,r)-l\\< [ \\AM{Ox,ey, s)\\ ds < f Qmax-i) ^ e'^^ r, (161) 

Jo v27r V 7 / 

where we used Lemma [3] to get the final inequality, which completes the proof. □ 

We now return to the proof of Lemma |8] 



25 



Proof of Lemma^ Using Lemma |9] and several triangle inequalities, we may approximate the unitary defined in (49 1: 

)Ul,{e,,ey + r, r)UY{0x + r, 9y,r) Ux{Ox-, dy, r) 

with the following version: 

V^^'\e^,ey,T) = 0^, r) C/j,(^,,)(0,, 0^ + r, r) UYiM){0. + r,0y, r) UxiM){0.,Oy,r) 

up to an exponentially small error. More formally, we have for all 9x,dy G [0,27r]: 



Vo{9,,,9y,r)-V^'''\9,,ey,r) 



< 256 



a7 
\/2^ 



7 y cr 



e 



(162) 



We introduce now the following useful unitaries that will aid us in the perturbative expansion of VQ^^\9x,9y,r): 

Fr{9x,9y, 81,82) — Uyf^j^j-^{9x,9y, 81) ul^^j^_j.^{9,^,9y + r, 82) Uy(m){Sx + 9y, si) Ux{M){0x,9y, 82) 
Gr{9x,9y, 81,82) = t^y(Af)(^2; + r, 9y, si) U\^^j^j-^{9x,9y + r, 82) Uy(m){9x + r, 9y, si) Ux{M){9x,9y + r, 82) 



where < 81,82 < r. Note that V, 



(M) 



(3 y'^xj'Jy 



,9y,r) = F, 



r) and that the following identities hold: 



9y, 81,82) — Fr{9x,9y,8i,0)Gr{9x,9y,8i,82) Fr{9x,9y,0,82) 



r\^X7 

Fr{9,,9y, 0,0) = Gr{9x,9y,0,82) = Gr{9x,9y,8i,0) ^ I 
dsGr{9x,9y,0,8) = dsGr{9x,9y, 8,0)^0. 

Finally, we define the following operators: 

D'^^'\9^,9y + 8) = Sj^'\H{9,,0,9y + 8,0),deH{9,,0,9,0)e=e,+s) 

r, 9y + si) = -i 



(163) 



^Y{9x.9y,8i,r) = iV^y^'\9,,9y + 8i)^iV^^'\ 



D 



X y'^xi'Jy 



+ 8) = Sj^'HH{9,,0,9y + 8,0),deH{9^,0,9,0)o=e,+s) 



Ax{9x,9y,82,r) = iV'^^'\9, + 82,9y)-iVf\ 



82, 9y + r) = -i 



+ 8',9y + Si)s' = s) d8 



h + 82, 9y + S')s' = s) d8 



Note that the support of Ax{9x,9y, 82,r) lies strictly within a 2M x 4M box, centered at the origin in the x — y 
orientation, whereas the support of Ay{9x, 9y, 82, r) lies strictly within a 4M x 2M box, centered at the origin. This 
is because each partial derivative eliminates terms that do not depend on the variables of interest and moreover, 
interaction terms composing the above operators have been trimmed to have a radius of support M. Our goal is 
to show that up to order 4 in r, the unitaries Fr{9^,9y, 8i,0), Gr{9x,9y, Si, 82) and Fr{9j.,9y,0, 82) are localized 
within the set flo, defined in (1191, and that the localized versions are rotations of the same operators evaluated at 
9x = 9y = 0. 

We begin by showing locahzation for Fr{9x,9y, si, 0). A similar argument applies to Fr{9x,9y,0, 82)- One can easily 
check that: 

dsFr{9x,9y, 8, 0)s=si = Ul^j^j^{9x,9y,8i) Ay{9x, 9y, si, r) Uy(m){9x, 9y, si) • Fr{9x,9y, si, 0) 

Noting that we only need to consider partials up to order 3 since \\AY{9x,9y, si, r)\\ = 0{r) (a discussion on how one 
may bound \\ds''D';^\9x + s', 9y + si)s'=s|| is given in the Appendix), we have: 

dsFr{9,,9y,8,0)s=o = AY{9^,9y,0,r) 



dlFr{9^^,9y,8,0)s=o = i[AY{9.„9y,0,r),D^^'\9,,9y) 



dsA 



d'^Fri9,,9y,8,0),=o = AY{9,,9y,0,r),D'/"{9,,9y) ,D'4">{9,,9y) 



's^Y\0x,0y, 



8,r)s=Q + Al-{9^,9y,0,r) 



+ ids 



AY{9x,9y, 8, r),D'^'\9„9y + s)! + d^s^Y{9x, 9y, 8, r),=o + 0{r^). 

J s=0 
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Torus T 



y=L 



y= U2 



[0.0) 




FIG. 3: Decomposition of the torus T into areas where the operators Sa^'^\H{O,6x,r,6y),d0H{6,Q,r,6y)\g^f)) and 
Sa^^\H{r,Ox,Q,Oy),deH{r,9x,9,Q)\e=o) act non-trivially. The dark shaded regions are defined as X = {s £ T : x{s) £ 
[-L/6,L/6] A y(s) G [-L/3 - 7?, L/3 + 7?]} and F = {s G T : y{s) G [-L/6,L/6] A a;(s) G [-L/3-i?, L/3 + -R1} and correspond 
to interaction terms ^c'^^^' (-^^(0, 0, r, 0), 6*, 0, r, 0)|e=o) with Z C X and 5c<*^>(-H"(r, 0, 0, 0), as<l>(Z; r, 0, 0)|e=o) with 

Z CY, respectively. On the other hand, the hght shaded regions X' = {s G T : x{s) G [-L/&, L/6] A y(s) G [2L/3, 4i/3]} and 
F' = {s G T : 2/(s) G [-L/6, L/6] A x{s) G [2L/3,4L/3]} correspond to terms 5^'"' (-^^(0, 0, 0, fy), ae$(Z; £», 0, 0, 6iy)|e=o) with 
^ C X' and ScS'*'\H{Q, 9^, 0, 0), dg^{Z; 0, Si, 0, 0)|e=o) with Z C K', respectively. Since each interaction has radius of support 
M < L/6, interactions in X' commute with interactions in Y and Y' and interactions in Y' further commute with those in X. 



Now, the crucial point is that the operators AyiOxjOy, s,r) and Dy^\6 
Arie^^ey, s, r),D'y^'\9x,0y + s)-Ry (Oy, Rx (e., i?^*'^(0, s) 



s) satisfy for M < L/2A: 



Ry [ey,Rx (e,, [Ay(0,0,s,r),Z?^*')(0,s) 



, D^^'\9x, 0y + s) - Ry (ey,Rx f^x, d\.^'\o, s) 



= 



which may be verified by considering the action of the twists Ryidy,-) and Rx{Ox,') on interaction terms with 
support contained in f2o- For example, note that for the truncated Hamiltonians defined in (71 1 with M < L/2, we 



have RYiOy,HMiZ;0,0,s,0)) = HM{Z;O,O,0y + s,0) for Z with support near ?/ = 1 and RxiO^, Hm{Z; s, 0,0,0)) = 
Hm{Z; 9x + s, 0, 0, 0)) for Z with support near x = 1. This follows from the trivial action of the truncated Hamiltonians 
Hm{Z; 0, 0, s, 0) and Hm{Z; s, 0, 0, 0) near the twists at ?/ = L/2 + 1 and x = L/2 + 1, respectively, and the fact that 
each Hamiltonian is the sum of interaction terms that commute with charge operators whose support covers the 
interaction. In general, the previous argument can be applied to any Hamiltonian that acts trivially on one of the two 
boundary lines of the charge operators Qx and/or Qy, effectively introducing a single boundary twist by applying 
the corresponding global twist Rx{9xi ■) and/or RyiOy, ■), on the whole Hamiltonian. 

Of course, in order for the previous argument to work, we need to show that each term we wish to rotate globally, 
has support away from the twists at x = L/2 + 1 and y — L/2 + I. To study the support of the terms derived 
in the above partials, we choose one of these terms as an illustrative example. In particular, we note that taking 
commutators of Ay (fl^,, s, r) with D^^\6,j.,0y + s), can only extend the support of AY{9x,0y, s,r) (a 4M x 2M 
box, centered at the origin) up to 2M in the y-direction, so terms affected by the twist on x = L/2 and y = L/2 
will not appear in the commutator for M < L/12, since then (2 + 2)M < L/2 — 2M. Finally, since the support of 

Dy^\6xj Sy + s) — Ry {Oy, Rx {^x, £'^^■'(0, s)^^ is away from the support of Ay{9x, dy, s, r), the commutator of the 
two operators is trivially 0. 

The above observations, combined with the partials we computed and the assumption < si, S2 < r, imply that 
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for a constant C > 0, the following bounds hold: 

\\Fr{9,,9y,si,0) ~ RYiey,Rx{0,,Fr{0,0,si,0)\\ < C {Q,,,^^aJ if r'> 
\\Fr{e^,ey,Q,S2) ~ RY{ey,Rx{e^,Fr{Q,Q,Q,S2)\\ < C {Q,,,^^aJ Lf r\ 



(164) 
(165) 



with the constant power of {QmsKOiJ L) coming from bounds on the norms of fifth-order commutators contai ning 
fQvvT^o Ti(^ha a si) and partials like ds'T^Y^^O^ + s' ,9y + Si)s'=s, whose norms may be bounded using (22 1 



terms like D 

and arguments similar to the one found in the Appendix, respectively. 

Moreover, the supports of the partial derivatives of F,.(0, 0, si, 0) and Fr{0, 0, 0, S2) up to order 4 in r, lie strictly 
within Hq, for M < L/48. By taking M to be a smaller fraction of L, we could have continued the expansion to 
higher orders in r, while keeping the supports of Fr{0, 0, si, 0) and Fr{0, 0, 0, S2) strictly within f^o- For our purposes, 
it suffices to consider errors of order 0{r^). 

We turn our focus, now, to showing the bound: 



5 „5 



\\Grie^,9y,s,,S2) - RYi0y,Rxi0., GriO, 0, Si, S2)))|| < C {Q^^^aJLf r 



(166) 



for some constant C > 0, while keeping the support of Gr(0, 0, si, S2) within f^o, up to order 4 in r (recalling that 
< si, S2 < r). We begin by setting: 



A 



Oy, Si, S2) = * C/i(M)(^x, Oy + r, S2) D\^'\e^ + r,ey + Si) Ux{M){Ox.ey + r, S2) - i D\^\e, + r, Oy + si) 



f^x(Af)(^- + [Df\e. + s, Oy + r),D^^'\e^ + r,0y + s,)\ Ux,^M){0.,Oy + r, s) ds 
A^x{Ox,0v,si,S2) = I Df\e,, + .S2,9y + r) - z + r, ey,si) Df\e, + S2,By + r) UYiM)iOx + r, By, si) 



C^F(M)(^- + ^' ^v^ [Dx'\dx + S2, &y + r),D\;'>{e., + r, 0^ + s)J {/y(M)(0, + r. By, s) ds 
Taking partials with respect to the variables si and S2, one may verify: 

ds'^Gr{9^,ey,s'^,S2)a'^=s^ = Gr{9x , By , Si, S2) ■ (^C/|,(j^^) (6*:^ + T, , Si ) Ay (6*3; , 6*2^ , Si , S2) UY(M)iBx+r,By,Si)j 

ds'^GriBx,By,si,s'2)s'^=s2 = (ul.^j^f^{Bx,By + r, S2) Axi^x, By, Si, S2) Ux(M){Bx,By + r,s2)^ ■ Gr{Bx,By,si,S2) 
Using the above relations, it is straightforward to check: 

dsiGr{Bx,By,Si,S2)si=Q — Ay (S^^ , 0y , 0, S2 ) 
ds^Gr{Bx,By, Si,S2)s2=0 = A.x{Bx,By, Si,0) 

dlGriB,, By, si, s2).,=o = i Wy{Bx, By, 0, s2),d'^'\Bx + r, By)] + (9,, A^(0„ 0^, si, sa)),^^^ + (A^(0„ 0^, 0, S2))' 



X, By, Si, S2)s2=0 



A 



X\"x 



, By,si,Q),Df\Bx,By + r) + (9,, A^(0,, si, S2)) + (A^(^?„ 0„ si, 0))^ 



We note here that we only need to evaluate the following partials at si = S2 = in order to have a complete picture 
of Gr{Bx, By, si, S2) up to order 4 in r: 



ds2idsiGr{Bx,By,Si,S2)si=o)s2=0 — ds^AyiBx, By, 0, S2)s2=0 



D 



ds2{dl^Gr{Bx,By,si,s2)si=a)^^^a = i 



D 



D 



(M) 
X 



(M) 



,9y + r), D^y^'^ (Bx + r, By)\ , D'^'' ' (B^ + r,By) 
{Bx, By + r), (^d,^D^y''\Bx + r. By + si)) 



^^'^'Bx,By + r),D^y^'\Bx + r,By) 



si=0 



Sl,S2) 



"2=0;,,^ 



si=0 



Df^ {Bx,9y + r), D'^f^ {Bx + r,By)\, of^ {B^ ,By + r) 
{d,,D^x\^- + S2.By + r)) , d'^^'\9x + r. By) 



dlA9lGr{Bx,By,si,S2)s2=o),. - ^ {dlA^x{Ox.By,s,,G)) ,Df\Bx,By + r) 



a,, (a2^A5,(0„0y,Si,O)) J +2 D'^'\Bx,By + r),D'4''\Bx + r,By) 

/ 6*2—0 L 
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Now, since [dg^Ax{0j:,9y, si,0)) = {d1^Gr{0x,0y,si,S2)si=o) we can see from Figure plthat all of the 
above commutators are rotated versions of the same commutators with 9x = 9y = 0, where we apply Ryidy, Rx{Ox, •)) 
to perform the rotation as in the study of Fr{0x,Oy, si,0). For example, we have for one of the fourth order terms 
implicit in the first term of the last line of partials evaluated above: 



D 



, 0y + r),D\^'\ex + r, ey)\ , d'^^'^Ox + r, 9y) 

where we have assumed that M < L/12 to make sure that all terms affected by the twists at a; = y = L/2 vanish. 
Furthermore, choosing M = L/48, guarantees that up to fourth order all commutators have support within fig- To 
see this, note that each new commutator includes terms with support at most 2M away from the support of each 
term in the commutator. For example, the above commutator has, potentially, the largest support of all the terms up 
to fourth or der. It is supported within a cross of radius 5M, with tw o axis of width 2M each, centered at the origin. 
Now, using ( |l62| with V^^\9x,9y,r) = Fri9 r,r), we get from ( |l63|l66| and several triangle inequalities: 



X 



.Df\9x,9y 



\\Vo{9x,9y,r) ~ RY{9y,Rx{9x.Vo{Q,Q.r)))\\ < C 



Q!7 

v/2^ 



7 



Moreover, setting W{r) to be the sum of terms up to order 4 in r in the Taylor expansion of the operator V^^ {0, 0, r) 
and noting that W{0) = 1 and drW{r)r=Q — 0, we have — 1|| < C {Q^a.->^otJ Lf' r"^ , for some constant C > 0, 

which completes the proof of this Lemma. □ 

Now, we come back to the proof of Lemma |6] which implies Corollary [5] 

Proof of Lemma^ Using Lemmas [t] and [s] with Ao,, — W{r) — 1, we get the bound ( 148 1 from the following estimate: 

I {^o,'i^o{ex,Oy,r)) - (vE'o,«'oM> I < I {^o,W{r)^o) - (^c^-oW) I 

+ |(^'o,M^(r)^'o) -(*o,^l^^[(0,0) ^ {9,x,9y)]RYi9y,Rx{9^,W{r)))V[iO,0) ^ {9^,9y)]^o)\ 
+ 1(^1/0,^^(0,0) ^ {9x,9y)]RY{9y,Rx{9xdW{r)~Vo{0,0,r)]))V[{0,0) ^ {9x,9y)]^o)\ 
+ |(v|/o,yn(0,0) -> {9x,9y)] [Vo{9x,9y,r) - RY{9y,Rx{9x,Vo{0,0,r)))] t/[(0,0) -> {9x,9y)]^o)\ 
< 2\\W{r) - T/o(0,0,r)|| + \\Vo{9x,9y,r) - RY{9y, Rx{9x,Vo{0,0,r)))\\ 

+ I {^0, {W{r) - 1)^0) - (*o, ^^^[(0, 0) ^ {9x,9y)] RY{9y, Rx{9x, {W{r) - 1))) V[{0, 0) ^ {9x,9y)]^o) \ , 
recalHng the bound || ^^(r) — 1|| < C {Qma^otJ L)^ r^, for some constant C > 0. □ 



PROOF OF MAIN THEOREM 



At this point, we can put everything together to prove the main result, the Quantization of the Hall Conductance: 

Proof of the Main Theorem. Going back to ( |128[ ) , we use the bounds derived in Corollary [s] Corollary |4] and Corol- 
lary [5] to complete the proof, noting that terms which are exponentially small in 0:^7^ dominate terms which are 
exponentially small in L. Finally, we set Gfij^j{L) — 0:7 with a given in the statement of LemmaW] □ 



DISCUSSION AND EXTENSIONS 



Hall conductance as an obstruction 



Given two different gapped Hamiltonians, Hq and Hi, it is natural to ask whether there is a path in parameter 
space {^^(s)}sG[o,i] which connects the two Hamiltonians such that H{s) obeys the same upper bounds R,J on the 
interaction range and strength, and given some lower bound on the spectral gap 7(5) of H{s), for ,s e [0,1]. As 
a corollary of our main theorem, we now show that there is an obstruction to finding such a path connecting two 
Hamiltonians with different Hall conductance, unless the gap becomes polynomially small in L, for some s S [0, 1]. 
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(0,0) (2jr,0) (0,0) (271,0) 

FIG. 4: Paths on the flux torus. Dashed hnes are used for the boundaries of the torus, and solid lines for paths, a) An example 
of two different paths that are compared in the first assumption. The first assumption states that the resulting projector is 
almost the same for both paths, b) Another example of two different paths, one of which travels vertically and then horizontally, 
while the other moves a slightly shorter distance vertically, then moves horizontally, then moves a short distance vertically to 
arrive at the same point. The second assumption describes a differential form of the assumption that these two paths give 
almost the same projector. 
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Suppose that Hq, Hi both have interaction range R and an energy ratio J/7 of order unity, with L » I. Then, the 
main theorem implies that the Hall conductance is very close to an integer multiple of e^//i in both cases. However, if 
the given integers, nQ,ni differ, then for any path, for some s G [0, 1], the Hall conductance must equal (uq ± 1/2)^. 

However, as a corollary of the main theorem, we have: 

Corollary 6. For any given R and Qmax, there exists a constant c such that for L > c( ln^^/^( J/7) 



< 



2h 



(167) 



Proof. For the given L, Gr.j,^{L) - {j/J)[{J/jy/Hn^^^^{J/j)/ln'^{J/-f)]^/^ - [In^/^^ j/^)]i/5 _ ini/2(j/^)^ ^hus, 
exp(— G|j j^/6) is inverse polynomial in L, and by appropriate choice of c the desired result follows. □ 

Thus, if riQ ^ ni, for some s, the gap 7(s) obeys 7(3) < c'JL"^/^ ■ ln"/^"(L), for some c' which depends on R, J. 



Fiber bundles 



We now relate our approach to Chern numbers. First, recall how the Chern number argument for the Hall conduc- 
tance worked. Under the assumption that the gap is non-vanishing everywhere, the Hamiltonian has a unique ground 
state density matrix \'i'o{dx,Sy)){'i'o{6x,dy)\ for all 9x,0y. This ground state density matrix is a rank-1 projection in 
a very high dimensional space (the Hilbert space of the whole system has dimension exponentially large in L^). 

Given this family of projections, they define a complex vector bundle, with the torus as the base space and C as 
the fiber. The connection given by adiabatic continuation is the canonical connection. 

One could also consider a different bundle, with the larger Hilbert space (that is, not just the ground state subspace, 
but all higher energy eigenspaces included) as the fiber. This bundle is trivial, being the direct product of the larger 
Hilbert space with the torus. In this case, if there is a unique ground state for all dx,Sy, then adiabatic continuation 
gives a connection that has the property of preserving the rank-1 projection onto the ground state. 

We now show that we can find a rank-1 projection that will be approximately preserved by quasi-adiabatic con- 
tinuation. We do this using general error terms (e below) so that the results can be applied to other quasi-adiabatic 
continuations. Our assumptions below correspond to assumptions that different paths starting at (0, 0) and ending at 
{(^xidy) give approximately the same state; the first assumption bounds the difference for two different paths, while 
the second bounds the change in the state for a slight change in the path, as shown in Figure |4] 

We obtain the desired projector by averaging over different paths; the reason this is necessary is that if we choose 
just one particular path from (0,0) to obtain the projector at (9x,0y), for example the path first increasing 9y and 
then increasing O^, this would give a discontinuity at 0^ = Oy = 27r. 
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Lemma 10. Suppose 'Dx{dx-0y),'I^Y{0x,6y) are quasi- adiabatic continuation operators and '^o is a state such that 

the following conditions hold: 

1. For any 9x,0y G [0,27r], let A", for a = 1,2,3,4 denote four different paths on the flux torus from (0,0) to 
{(^xj(^y) follows. We can move from (0, 0) to (0, dy) by either increasing Oy from zero or by decreasing it from 
27r, and, then we move from {0,0y) to {6x,0y) again by either increasing Ox from, zero or by decreasing it from, 
2iT. The paths with a = 1, 2 correspond to the paths in which we increase Oy, while those with o = 3,4 correspond 
to those with decreasing Oy, while the paths with a = 1,3 correspond to those with increasing Ox and those with 
a = 2,4 correspond to those with decreasing Ox- Let Ua, for a = 1,2,3,4 denote the unitary evolution using 
quasi- adiabatic continuation along each of the corresponding paths. Then, we require that 

|(*o,C/if/6*o)|'>l-e', (168) 

for all a, b, for sufficiently small e. 

2. Next, let UyiOxTOy) denote the unitary evolution corresponding to the vertical segment of path [/" (the evolution 
in Oy), while UxiOx^Oy) denotes the unitary evolution corresponding to the horizontal segment (the evolution in 
Ox). Note that Uy depends only on Oy, not on Ox. Then, we require that 

\\i{VYiOx,Oy)U^{0x,9y) - U^{0x,0y)VY{0,0y)U^{0x,ey)} \^o}{^o\U''{9x,0y) + h.c.\\ < e. (169) 

3. The operators VxiOx, 0y),'Dx{0x,0y) are differentiate with respect to Ox,Oy everywhere and are periodic in 
Ox,Oy with period 2tt. The norms WDxiOxjOy)]], \\T>Y{0x,0y)\\ are bounded by some constant D. 

Then, there exists a rank-1 projector, P{9x,0y) which is periodic in 9x,0y with period 27r such that 

||P(0,0)-|*o)(*o|||<Ce, (170) 

and such that 

XI 9y), P{Ox, 0j/)] — dg^P{Ox, 9y)\\ ^ C(l + D)e. (171) 

and 

\\i[DY{9x,ey),P{Ox,ey)]-de^P{0x,ey)\\ < C{l + D)e. (172) 
Proof For a = 1,2,3,4 let P{9x,9yY = Uo\^o){^o\Ul. By the first assumption, 

Tr(P(e,, eyrP{9x, Oyf) >l-e\ (173) 

For a e {1,2,3,4}, and Ox,Oy e [0,27r], let w°-{0x,0y) = w^{0x)wy{0y) denote the weight of the path A^, where 
Wy{Oy) is equal to 0y/2iT if the path initially moves by decreasing Oy from 27r and is equal to 1 — 9y/2'K if the path 
instead initially moves by increasing 9y from and w'^{9x) is equal to 9x12-^ if the path then moves by decreasing 9x 
from 2-K and is equal to 1 — 9x12^ if the path instead moves by increasing Ox from 0. 

For Ox,Oy G [0,27r], let 0{Ox,Oy) - ELi ^«"(^- ^y)^(^- ^j/)"- 

Since 0{2tt, 9y) = 0(0, 9y) and 0{9x,Q) = 0{9x, 27r), we can extend the definition of 0{9x, 9y) to a periodic function 
for all 9x,9y. The operator O is differentiable everywhere except when 9x or 9y is an integer multiple of 2n. In the 
rest of the proof, whenever we evaluate a partial derivative of 0{9x, 9y), our bounds are valid for all 9x, 9y except for 
the lines where 9x or 9y is an integer multiple of 27r. 

First, we bound the norm of i['Dx{9x, 9y),0{9x, 9y)] — dg^O{9x, 9y). Note that for each a, 

i[Vx{9x, 9y),P{9x, 9yr] = de^P{9x, 9yr. 

Thus, 

4 

i[Vx{9x,9y),0{9x,9y)]-deMSco,9y) = -J2ide.^"i0x,9y))P{9x,9yr = 

o=l 

{t " (^(^-^^)' - Pi^-^^y)') + ^ {m,Oyr - P{0x,9yf) . 
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Using II All < II All 2 = v/Tr(A2), true for any Hermitian operator A, we have: 

¥['Dx{o.,0y),o{0x,ey)] - deM0.,Oy)\\ < \\P{o.,Oyy - p{ex,ey)^ + ^\\p{ex,9y)^ - p{ex,ey)' 



^ max {||P(0,, - P{0x, 0y)^\\^ , \\P{ex,Oyf - P{9x, 0y)%} < |^ e, 



where we used (1731 to get max | ||P(6I^, 6ly)i - P{0x,0y)^\\l ,\\P{0x,0yf - Pi0x,Ov)*\\l^ < 21^, for the last bound. 

We next evaluate i['DY{0x,()y),O{0x,0y)] — dg^O{0x,0y)- Now the difference between quasi-adiabatic con- 
tinuation and partial derivative, namely i['DY{0x,0y), P{0x,0y)°'] — dg P{0x,0y)'^, is non- vanishing. The op- 
erator norm of this difference is equal to \\ {iVY{0x,0y)U''{ex,9y) - ds^U%0x,Oy)} \'^o){'^o\U''{0x,0y) + h.c.\\. 
This is equal to \\l{VY{0x.0y)U^{dx,0y) -U^i0x,ey)VYiO,0y)U^{0x,0y)}\'fo){\U''{9x,0y) + h.c.\\ < CI. So, 
\\i[VY{0x.0y).O{0x,0y)] - de^O{ex,0y)\\ < Ce+ {l/27r)ma.x{\\P{0x,0y)^ - P{0x,0yn,\\P{0x,0yf - P(0.,0,f ||}, 
which is bounded by Ce. 

Given the operator O(0x, 0y), we can define a smoothed operator by convolving it: 

6{0x, 9y)= J J O(0x, <l)y)H(j)x - 0x, <i>y " 0y)d(bxd(j)y , (174) 

where k{x,y) has support in a small square near x — y — Q, k is positive, and the integral of k is equal to unity. 
For differentiable by choosing the size of the support of k sufficiently small, we can construct an operator 

0{9x,9y) that is everywhere differentiable, is periodic, and obeys the bounds 

\\i[Vxiex,0y),diex,0y)]- deMOx,0y)\\ <Ce, (175) 
\\i[VY{0x,Oy),di0x,0y)]- de^d{0x,0y)\\ < Ce, (176) 
\\d{9x,0y)-O{0x,0y)\\<Ce, (177) 

for all 0x,0y 

We now define P{0x,0y) to be the projector onto the eigenvector of 0{6x,dy) with the largest eigenvalue. By 
( [l73| , we have Tr{d{0x,0y)P{0x,9yT) > 1 - Ce for a e {1,2,3,4}. Hence, by definition the eigenvector P{0x,0y) 
has eigenvalue greater than or equal to 1 — Ce. This eigenvector is unique for sufficiently small e and the next 
largest eigenvalue is less than or equal to Ce. We now relate these bounds to bounds on \\i[DY{0x,9y), P{Sx,9y)] — 
dg^P{0x,0y)\\- For notational convenience, we write O for O{0x,0y) and P for P{0x,0y) when no derivative is taken. 
Since ||P — 0|| < Ce, we have 

\\[VYi0x, 0y),P{ex,0y)] - [VY{,ex,Oy),d{9x..0y)]\\ < 2\\Vy{0x, 0y)\\Cl < DCe. 

Let X{0x,0y) be the largest eigenvalue of O{0x,dy) and v{9x,0y) be the corresponding eigenvector so that P = 
X{0x,0y)-^H0x,0y)){vi0x,9y)\- Then, we have 

(1 - P)deX9^,0y) - \ - (1 - P) (de^di0x, 0y)) v{0x,0y), (178) 

Then, since all eigenvalues of O other than A are close to zero, we have that || (A(6'a;, 0j,) — 0)^^(1 — P) — (1 — P)|| < Ce. 
Thus, 

\\deP{0x.0y) - X{0x.0yY^ {(1 - P)dep{0x,0y)P + Pdep{0x,0y){\~P)\ II < Ce. 

Since dop{0XT0y) is close to i\VY{0x,0y),O{0XT0y)] which in turn is close to i\VY{0XT0y)TP{0x,0y)], and since A~^ 
is close to unity, we find that 

\\de,P{0x, 0y) - « {(1 - P)[VY{0x,Oy),P]P + P[Vy {9 x , y) , P]{1 - P)} II < C(l + D)e. (179) 

However, (1 - P) [I?y(^^„ 0^), P]P + P[I?y(0„ 0„), P](l - P) = [Py ((?„(?,), P] so 

\\de^P{0x,0y) - i[VY{0x,0y),P]\\ <C{\ + D)e. (180) 

The proof for the bound on \\i\Dx{0x, 0y), P{9x, 0y)] — de^P{0x, 0y)\\ is similar. □ 
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Having defined the projector P{9x, dy), we can define a fiber bundle with C as the fiber, and the projectors define 
a canonical connection. One can show that the Chern number with this connection is close to the Hall conductance 
in units of e^//i. 

Fractional Hall Effect 

We now discuss the fractional Hall effect. We do not claim to present a complete proof of the fractional quantization 
at the same level of detail as in the integer case; rather, our goal is to sketch a proof to indicate where modifications 
need to be made in comparison to the integer case. We will present a detailed proof elsewhere. In the case of the 
fractional quantum Hall effect, we consider the same class of Hamiltonians Hq, but now we assume that there are 
q different eigenstates of Hq, for a — l---q, which we refer to as ground states. We assume that the energy gap 
between the highest energy ground state and the rest of the spectrum is at least 7. 

The fractional quantum Hall effect can also be treated using our approach, with one additional assumption about 
topological quantum order required and with an appropriate definition of the Hall conductance, as discussed below. 
We first need to define topological quantum order. Following [IT , we define that the q ground states have topological 
quantum order with error if, for any operator Oioc with ||0(oc|| — 1 supported on a set with diameter I or less, 
we have that the q-hy-q matrix M with matrix elements 

A^b- («'S,0/oc«'o), (181) 

obeys \\Mab — cSab\\ < e for some constant c. That is, the projection of Oioc into the ground state sector must give 
an operator close in operator norm to a multiple of the identity operator. Given this definition, we require that the 
system has {L/2 — l,e) topological order for suflicicntly small e. 

Next, we need to correctly define the Hall conductance. This problem was physically discussed by Thouless and 
Gefen[r5]. They pointed out that, at experimentally relevant temperatures and voltages, the terms in the Kubo 
formula involving matrix elements between the different ground states should be neglected. Thus, we should define 
the Hall conductance in any normalized state ^Pq which is a linear combination of ground states by 

a.,(vI/o) = 2Im|^J-vl/o(O,0,),„=o,(l-P,)J^*o(^.,O)e,=o^|- (^27r^^, (182) 

where Pg is the projector onto the ground state subspace: Pg — J2l=i l^o)(^ol- We will see below that quasi-adiabatic 
continuation, with an appropriate choice of a, provides a mathematical realization of the idea of Thouless and Gefen 
that, in the fractional Hall effect, experimental measurements are performed slowly enough to avoid exciting states 
above some large energy gap, but quickly enough that we "shoot through" any avoided crossing of the low energy 
states [131 [TB]. 

We will see that this definition gives a Hall conductance which is almost independent of the particular "ifo we choose. 

We will define axy to be the average of axyi'^o) over a = 1, q. Using these assumptions, we will prove quantization, 

2 

up to small error, of the Hall conductance in units of with q the dimension of the ground state subspace. This 
topological order assumption is certainly necessary to prove this. To see why the topological order assumption is 
necessary, let us consider an example of a system lacking this assumption. Imagine that a system is close to a first 
order phase transition between a state with filling fraction v — 1/3 and another one with ly = 1/5. By tuning the 
parameters in the Hamiltonian, we can construct a situation in which the two phases have the same energy and the 
ground state is 8-fold degenerate. However, in 3 of the ground states, we will observe a conductance close to 1/3 and 
in 5 of them it will be close to 1/5, neither of which is a multiple of 1/8. 

We will use the same choices of a,r as in the integer case. However, throughout this subsection we use o(l) to 
denote quantities which vanish at fixed i?, J, 7 in the limit of large L and vanishing e, while we use 0(e) to denote 
quantities which are bounded by a constant times e in the limit of small e at fixed R, J, L, 7, q, qmax- Such quantities 
which are 0(e) will depend on R, J/7, L, q, qmax only polynomially, and hence if e is exponentially small in L will be 
small. Due to the detail involved in the estimates in the integer case, here we will be content to sketch the errors, 
although, for sufficiently small e, the final error bounds on quantization (now up to a multiple of /hq) will be the 
same. We use "close" to denote that two quantities are within o(l) of each other. 

The assumption of topological order implies that the energy difference between the ground states is at most eL^J, 
since the Hamiltonian is a sum of local operators, each with energy difference bounded by eJ. Given this energy 
difference, it immediately follows [T^ that the matrix elements of the quasi-adiabatic evolution operators 'Dx,'Dy 
between the different ground states a,t 9x = Oy = are also 0{e). An alternate way to bound the matrix elements of 
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Vx:T^Y between the different ground states is to note that each one is a sum of operators Sa{H, Az) and each such 
operator, up to corrections which are exponentially small in L, is local up to length scale L/2 — 1. 

One can now repeat the integer proof, with one important difference described below. First, one can prove that 
evolution of a state ^'o which is a linear combination of ground states around the small loop of size r near the 
origin gives a phase which approximates e^^^""*^*"-'^'' There are two sources of error here. The first is 

that quasi-adiabatic continuation produces non-vanishing matrix elements between the ground states, while we have 
defined the Hall conductance with the projector (1 — Pg) which removes those intermediate states. However, quasi- 
adiabatic continuation has only small matrix elements between the ground states and these terms are 0(e) for any 
fixed R,J,L,j,q,qmax which in turn implies a fixed value of a. The second source of error is the same as in the 
integer case: errors due to the non- vanishing value of r, and errors due to the finite value of a meaning that we only 
approximate the quantity (1 — Pg)de^'^ o{0xiO) by quasi-adiabatic continuation. 

Since the evolution around the small loop can be approximated by a local operator, it is close to a constant 
times the identity operator when projected into the ground state subspace, and hence that (Jxyi'^o) is approximately 
independent of ^'q, as claimed. 

The proof that the product of evolutions around the small loops is equal to the evolution around the big loop is 
identical to the previous case. 

The energy estimates and the translation lemma both rely on the partial trace estimates. The partial trace estimates 
can be repeated as before, with the same error estimates as in the integer case, plus errors which are 0{e). In proving 

(2) 

the partial trace estimates in the fractional case, note that the product exp{—iQxd)Ux (0,0,9) is close to a phase 
multiplied by the identity operator when projected into the ground state subspace. This can be shown by using the 
topological order assumption: since the Hamiltonians H(9, —6, 0, 0) are unitarily equivalent for all 9, we need to show 
that the partial derivative of exp{—iQx9)Uj^ (0,0,9) is close to a multiple of identity operator after projection into 
the ground state subspace; however, this follows because this partial derivative is a sum of local operators at 9 = 0. 

The translation lemma [7] can be repeated as before, with the same error estimates as in the integer case, plus errors 
which are 0(e). Thus, we can, with the same error estimates up to terms which are 0{e), show that the analogue of 
corollary [5] holds, where now ^'q in corollary [5] can be replaced with any linear combination of the ^'q, a = l,...,q, 
and ^o('') is defined by the quasi-adiabatic continuation of that state: 

W) = ^0(0, 0, r) |vl/o(0, 0)) = 4(0, 0, r) [7^(0, r, r) Uyir, 0, r) Ux{0, 0, r) |vl/g(0, 0)) . (183) 

The evolution around the big loop can again be exactly decomposed into a product of evolutions F^[(0, 0] — ^ 
{9x,9y)]V0{9x,9y,r)V[{O,O) {9x,9y)] around small loops of size r. Since the evolution around each small loop is 
local and so is close to e*'^"'"'^'' /^^^('^ times the identity operator in the ground state subspace, this implies that 
the unitary V(3(0,0,27r) is close to a multiple of identity in the ground state subspace. That is, Zab = (^O' ^o(^'''')) 
satisfies \\Z — zl\\ < o(l), for some phase z, with z close to g^'^^y'^'^h/e^ ^ 

However, here is the difference. Unlike the integer case, we do not know that Z is close to the identity matrix 
itself, for the following reason. The energy estimates given in Proposition [l] for evolution around one side of the 
large loop can again be repeated as before, with the same error estimates up to terms which are 0{e), to show that 
quasi-adiabatic evolution of a ground state around that side gives a state which is close to a linear combination of 
ground states. However in this case the operator Ux{0,0,2'k) is supported on a set of diameter L, so the assumption 
of topological quantum order does not tell us that this operator is close to constant times the identity operator in the 
ground state sector. 

Since we no longer know that the phase cancels between evolution around opposite sides of the big loop, we must 
determine what this phase is. Let and Uy denote q x q matrices defined as follows: 

{u.x)ab^{%,UxiO,0,27r)^'„), iuy)ab^{%,UYiO,0,27T)^Sj'^,). (184) 

Then, Z — v)yu\. Uy Ux plus error terms which are due to leakage out of the ground state sector. In particular, the 
product of operators u\, Uy corresponds to projecting back into the ground state sector after each side of the 
big loop, hence the error terms which come from projecting off the ground state are small due to Proposition [T] in a 
manner similar to the proof of Corrolary|4| Since zl|| < oil) for some phase z, we have \\v\j u\. Uy Ux — zl\\ < o(l). 
Moreover, we have that Hu^^Ma; — 1|| < o(l), — 1|| < o(l) since, for example, 

{ulux)ab = (^o, Uj^iO, 0, 2TT)PoUxiO, 0, 2^X) - Sab - C/i(0, 0, 2^)(1 - Po)C/x(0, 0, 27r)^'o) , 

with Pq the projection onto the ground state sector and the second term above is small due to the energy estimate 
in Proposition 111 applied to the quasi-adiabatic continuation with Ux{0,0,2'k) of ^'q and ^f^j, separately. Finally, we 
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note that for q x q matrices A,B satisfying B = zl + A, with ||^|| < o(l), we have |det(B) — z'^j < o(l). Putting 
everything together and using the triangle inequality, we get from det(uj^ u^ Uy Ux) = det(u^Ux) ■ det(uj^Uy) and 
max{|det(u1,Ux) - 1|, |det(u1;uy) - 1|} < o(l): 



|det(uj^ u^ Uy Ux) — z'^l < 0(1) 



|det(u^ Ux ) •det(u|;uy) -z^l < 0(1) 



|l-z^|<o(l). 



Hence, z is close to a q-ih root of unity and hence axy is close to an integer multiple of {l/q)e^ /h. 
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APPENDIX 



ESTIMATING THE DIFFERENCE OF POWERS 



In this section, we prove the implication ( 131 1. First, we have 6 > |e | — |6 — e | > 1 — e> 1/2 and using induction 
it is easy to show that 6™ > 1 ~ me. Moreover, using Euler's formula e'^ = cosO + ismO, we have that 

|&- e'^P = (1 - 6)^ + 26(1 - cos 61) = (6 - cos e)^ + sin^ 6I < =^ sm\e\ < e, cos|6i| > y/l^ > ^3/2, 



which implies that cos^ is decreasing for (f> e [0, |0|] and hence, sin 1 6*1 = /d^' cos 



sin|6l| > 16*1 cos |6'| 



\e\ < 2e/V^- Moreover, 1 - cos(to6') = /g"^ # /^f cos 77^77 < /"^ |(/)| d0 < \e\^/2 < fm^e^. Armed with these 
bounds, we get the inequaUty |6™ — e™®| < W| me as follows: 



9 7 



= (1 - + 26™ (1 - cos(m6')) < - m^ t\ 

3 



and taking square roots completes the proof. 



Bounding the norm of higher partials 

The following general formula for a differentiablc family of Hamiltonians H(6), which can be verified by differenti- 
ating both sides with respect to t, will be instrumental in the analysis that follows: 

(185) 



\ J 4>=9 \ 1 4>=e Jq 

Using the above identities, we derive: 

{r„«(*'°'^-°)(5«ff (0, 0, 9, 0))o^o, }^^^^ - r^^io.fi.o^.o) (^(^g^ (5^ ^(^^ q, 6, 0)),^,^ 
|(9^e™^(^^°-«-"))^_^ e-™^(«^-"^^-")| ,r„^(^^^O-«-")(ae7?(0„O,0,O)e=eJ , 

(9^e™^(^-"^^'"))^_^ e-™^(^-"^^-")|,r„^(^-°'«-OH99i?(0.,O,0,O)e=eJ 
From ( |185| we have: 

V / (t>=6m Jo 

V / <f>=6,, ./n 



(186) 



(187) 
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Using the above formulas with the definition of T>y{0x, 0y), we get: 



< 



\t\Sa{t)dt 



{d^VYie^A)), 



< 



{de H{e,, 0, e, 0)),^,^ {d^ H{ct>, O, Oy, 0))^^, 



(9277(0,, 0,0,0)) 



Sa{t)dt 



t\Sa,it)dt]+ (9^i/(0„O,0,O))^ 



f Sa{t)dt 



with similar bounds for partials involving l?x(^x,0?/)- Now, recalling that Sa(t) is a Gaussian, we have 
J^^\t\ Sa{t) dt = and J^^t^Sa{t)dt = c? . Finally, using estimates similar to the ones derived in (SjjOl, it 
should be clear that the above partials have norms bounded by G {QmaxCtJ L)"^, for some C > 0. In general, one could 
show that TO-th order partials of VYiOxjdy) and 'DxiOxj&y) would have norms bounded by C {Qyna^ctJ L)™^^ . The 
same argument applies to partials of ^^^■'(fl^,, 6y) and 'D'^^^\9x, 9y), after we decompose each into a sum of individual 
interactions. 



DECOMPOSING THE EVOLUTION, ONE LOOP AT A TIME. 

Let N = {2n)/r and define C/(jv-m)+nW — V^{m ■ r,n ■ r)V(^{m ■ r,n ■ r,r)V{in ■ r,n ■ r), for m,n g [0, — 1] to 
be the evolution operator corresponding to the cyclic path leading to and around the square with lower-left corner at 
(m ■ r,n ■ r) in flux-space (see Fig.|2|. Then, we have the following identity: 



Uo(0, 0, 2^) = Un2 [/jv2_i ■■■U3U2U1 



(188) 



as can be verified by the decomposition described in Fig. [2] To help us track the contribution from each term in the 
decomposition as we insert Pq — |^'o)(^ol and Qo = 1 ^ -Po after each cyclic evolution, we introduce the scalars 



P[s,r] (*0,C/rC/r-l---C/s*o), 



and 



{^0\ UrUr-l---Us-lQoUs\^o) 



Moreover, we note that pr = (^'o,?7r^'o) and in particular, pN = (^'o, ^'o(0i Oi ^)) = (^07^'o(^))- The quantity we 
are trying to bound in ( 148 I is 



(vl/o,vl/o(27r))-(vl/o,vl/o(r))(^)^ 
Using the triangle inequality, we get: 



P[im'^]-{pn) <\p[i,m]-PiP2---PN^-iPN^\+ PiP2---Pn^-iPn^ - (pn) 



(189) 



so it suffices to bound each term on the right. We start with a bound for |p[i.jv2] ~ PiP2' ' ' Pn^-iPn^\- Using the 
following recursive relationship: 

P[l,Ar2] = PlP[2,N^] + Q[1,N^] 
P[2,N2] = P2P[3,N^] + q[2,N^] 



we may write: 



P[i,m]- PiP2 - ■ ■pm-iPN'^ ^ X! ^'i ■ ■ 



(190) 



We turn our attention to jgjj jy^]!- We begin by observing that the following series of bounds hold: 

\q[^,N-]\ = I (^-0, C/JV- • • C/.+ lQoC/»*o) I < WQoU^ 1*0) II - - \P^\^ < V'^i^ - \P^\) < ^2(1 - \Pn\) ' 



Pn\ 
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Since, we have iV^ — 1 terms in ( 190 1, the triangle inequaHty and the above bound together with the fact that \pi\ < 1 
for all i > 1, imply: 



\P[l,N^] -PlP2 ■ ■ -Pn^-iPn^I < {N'^ - 1) \/2(l - \pn\) + sup y^2\pi-pN\ 

Moreover, assuming supjgj;^ jy^] \Pi — Pn\ = S/N'^ for some S € [0, 1], we get: 



(191) 



(pn)^^ -piP2---PN^-iPN^ <{l + S/N^)^ -l<e'^-l< / eydy<e'^6<e-d. 

Jo 

where the first inequality follows from expanding the product 

P1P2 ■ ■■pm-ipm = [PN + [pi -Pn]) ■■■{pn + [pm-i -Pn]) [pn + [pm -Pn]) 

and then, after subtracting the term (pn)^ , using the triangle inequality along with \pn\ < 1. The second inequality 
follows from comparing the binomial expansion of (1 + x/m)"^ term by term with the Taylor expansion of and 
using the simple inequality (™) < Putting everything together, we get the final bound: 



(192) 



('i'o,«'o(2^))-(*o,*oM)^^)^ 



<N^ ^2{1 - \pn\) + J2 sup \p^-pN\ + e sup \p^-pN\] (193) 



recalling that N = ^ and noting that 1 — |pAr| is bounded in (1291, while supjgjj^ ^2] \pi — Pn\ is bounded in (1521 
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